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Section  1 


IN1R0DUC110N 


1 . 1  Pr^  Jim  Pjurgo^e 

Computer  program  HMR5?  computes  basin-average  precipitation  for  Probable 
Maximum  Storms  (PHS)  In  accordance  with  the  criteria  specified  In 
Hydrometeorological  Report  No.  S2  (National  Heather  Service,  198?).  That 
Hydrometeorological  Report  (HMR)  describes  a  procedure  for  developing  a 
temporal  and  spatial  storm  pattern  to  be  associated  with  the  Probable  Maximum 
Precipitation  (PMP)  estimates  provided  In  Hydrometeorological  Report  No.  51, 
"Probable  Maximum  Precipitation  Estimates  -  United  States  East  of  the  105th 
Meridian."  The  U.S.  National  Heather  Service  (NHS)  has  determined  the 
application  criteria  1n  a  cooperative  effort  with  the  U.S.  Array  Corps  of 
Engineers  and  the  U.S.  Bureau  of  Reclamation. 

Other  reports,  HMR  Nos.  36,  43,  49  and  55  (NHS,  1961,  1966,  1977,  and 
1983,  respectively)  describe  the  PMP  In  other  regions  of  the  U.S.,  Fig  1. 

This  program,  HMR5?,  Is  applicable  only  to  the  eastern  U.S.,  and  Is  Intended 
for  areas  of  10  to  20,000  m1».  (HMR  Mo.  52  also  contains  a  l-ro1^, 

1-hr  PMP).  A  time  Interval  as  small  as  5  minutes  can  be  used  for  storm 
definition.  Before  using  the  HMR5?  program,  one  should  be  thoroughly  familiar 
with  the  procedures  described  In  Hydrometeorological  Report  No.  52. 


figure  1.  Regions  Covered  by  Generalised  PMP  Studies  (NHS,  1980) 
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The  generalized  PMP  maps  of  HMR  No.  51  are  stippled  In  two  regions. 
Indicating  estimates  may  be  deficient  because  of  orographic  Influences.  Hajor 
projects  within  the  stippled  area  should  be  considered  on  a  case-by-case  basis 
and  expert  hydrometeorological  guidance  should  be  sought. 

Data  required  for  application  of  the  HHR5?  program  are. 

•  X.Y  coordinates  describing  the  river  basin 

and  subbasin  watershed  boundaries; 

>  PMP  from  HMR  No.  51  (NWS.  1978);  and 

•  Storm  orientation,  size,  centering,  and  timing. 

Ihe  program  computes  the  spatially  averaged  PMP  for  any  of  the  subbasins  or 
combinations  thereof.  The  Probable  Maximum  Flood  (PMf)  can  then  be  computed 
as  the  runoff  from  the  PMS,  using  an  appropriate  precipitation-runoff  program 
such  as  HEC-1  (Hydrologic  Engineering  Center.  HEC,  1981).  A  typical 
application  of  HMR5?  does  not  produce  a  PMS.  The  PMS  Is  defined  by  the  Corps 
of  Engineers  to  be  that  storm  which  produces  the  PME .  Thus,  the  PMS  can  only 
be  determined  by  computing  (and  maximizing)  runoff.  That  Is.  the  runoff 
characteristics  of  a  watershed  must  be  considered  In  PMf  (and  therefore  PMS) 
development . 

HMR  No.  52  requires  that  a  critical  storm-area  size,  orientation, 
centering  and  timing  be  determined  which  produces  the  maximum  precipitation. 
Ihe  HMR52  computer  program  will  optimize  the  storm-area  size  and  orientation 
In  order  to  produce  the  maximum  basin-average  precipitation.  The  user  must 
provide  the  desired  centering  although  the  centroid  of  the  basin  area  Is 
provided  as  a  default  option. 

Ihe  user  must  specify  the  time  distribution  for  that  storm.  Using  that 
time  distribution  Information,  the  HMR52  program  will  produce  a  data  file 
containing  the  Incremental  basin-average  precipitation  values  for  every 
subbasin  requested.  That  precipitation  data  file  will  subsequently  be  Input 
to  a  rainfal  1 -runoff  model,  such  as  HEC-1,  for  computation  of  the  resulting 
flood.  The  user  then  analyzes  the  various  storm  variables  and  recomputes  the 
floods  In  order  to  determine  the  storm  which  produces  the  maximum  runoff. 

That  storm  and  runoff  are  defined  as  the  PMS  and  PMF,  respectively. 

1 . 2  Computer  Requirements 

The  HMR52  computer  program  requires  a  computer  with  45K  (decimal)  words 
of  core  storage  and  7  scratch  tape/disk  files.  Plots  of  the  basin  geometry 
and  storm  patterns  can  be  made  on  a  line  printer.  Section  10  of  this  manual 
specifies  detailed  computer  hardware  and  software  requirements. 

1 . 3  Acknowledgements 

Ihe  computer  program  HMR52  was  written  by  Paul  B.  Ely  of  the  HEC. 

John  C.  Peters  provided  much  valuable  assistance  in  the  design  of  the 
program's  capabilities  and  applications  methodology. 


Section  2 


PROBABLF  MAXIMUM  S10RM  ANALYSIS  PROCEDURE 


? • 1  Probable  Maximum  Precipitation  Definition 

Probable  Maximum  Precipitation  (PMP)  Is  theoretically  the  greatest  depth 
of  precipitation  for  a  given  duration  that  Is  physically  possible  over  a  given 
sUe  storm  area  at  a  particular  geographical  location  at  a  certain  time  of  the 
year.  Hydrometeorological  Report  No.  51  (HMR  No.  51)  contains  generalized 
(for  any  storm  area)  all-season  estimates  of  PMP  for  the  United  States,  east 
of  the  105th  meridian.  Fig.  1. 

? . 2  Probable  Maxi mum  Storm  Definition 

Probable  Maximum  Storm  (PMS)  Is  a  hypothetical  storm  which  produces  the 
Probable  Maximum  Flood  from  a  particular  drainage  basin.  Hydrometeorological 
Report  No.  52  (HMR  No.  52)  provides  criteria  and  step-by-step  Instructions  for 
configuring  a  PMS  using  PMP  estimates  from  HMR  No.  51.  Key  concepts  upon 
which  the  procedures  In  HMR  No.  52  are  based  are  as  follows. 

2.2.1  S£a tial  Distribution 

Ihe  spatial  distribution  of  the  PMP  Is  governed  by  principals  described 
under  four  headings:  Isohyetal  shape,  orientation,  storm-area  size,  and 
spatial  variability. 

(U  Isohyetal  shape.  The  PMS  Is  represented  by  elliptical  Isohyets, 
each  of  which  has  a  ratio  of  major  axis  to  minor  axis  of  2.5  to  1.  Standard 
ellipses  have  been  established  containing  areas  from  ID  to  60, ODD  ml* 


Figure  2.  Standard  Isohyetal  Pattern  (NWS,  198?) 


(II)  Orientation.  There  Is  a  preferred  orientation  for  storms  at  a 
particular  geographic  location.  That  orientation  ls  related  to  the  general 
movement  of  storm  systems  and  the  direction  of  moisture-bearing  winds. 

Contours  of  preferred  orientation  are  shown  In  Fig.  3.  When  developing  a  PMS, 
It  Is  generally  desirable  to  orient  the  storm  to  produce  maximum  precipitation 
volume  In  the  watershed.  PHP  will  be  reduced  by  an  adjustment  factor  (shown 
In  Fig.  4)  when  the  storm  orientation  differs  from  the  preferred  orientation 
by  more  than  *  40  degrees. 

(III)  Storm-Area  Size.  The  average  precipitation  depth  over  an  area  1s 
PMP  for  one  and  only  one  area  size.  This  Is  the  "storm-area  size."  The 
average  precipitation  on  areas  larger  or  smaller  than  the  storm-area  size  Is 
less  than  PHP  for  the  larger  or  smaller  areas.  Fig.  5  Illustrates  this 
concept  for  a  storm  area  of  1000  ml^.  The  storm-area  size  Is  chosen  to 
yield  the  maximum  precipitation  volume  from  a  given  drainage  basin. 


Figure  3.  Preferred  Orientation  for  PMS  (NWS,  198?) 
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Figure  4.  PMP  Orientation  Adjustment  Factors  {NWS,  198?) 
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Figure  5.  Comparison  of  PHP  Depth-Area  Relation 
with  1,000  ra1»  PMS  (NWS,  198?) 

(1v)  Spatial  Variability.  Spatial  variation  of  precipitation  Is  a 
maximum  during  the  6-hr  period  when  the  maximum  precipitation  occurs.  Spatial 
variation  diminishes  for  the  second  and  third  largest  6-hr  amounts.  For  the 
remaining  6-hr  periods,  the  wlthln-storm  precipitation  Is  uniform,  but  there 
Is  spatial  variation  In  the  residual  precipitation  occurring  outside  the 
elliptical  boundary  that  corresponds  to  the  storm-area  size.  HMR  No.  52 
contains  nomograms  which  express  spatial  variation  for  each  6-hr  period  as  a 
percent  of  PMP.  Percentages  for  selected  area  sizes  are  tabulated  In  Tables  1 
through  4.  For  each  Isohyet,  the  percent  of  PHP  for  the  storm  area  Is 
Interpolated  from  Tables  1-4.  Ihose  percentages  are  multiplied  times  the  PHP 
to  obtain  precipitation  for  each  Isohyet  for  each  6-hr  Interval.  A  graphical 
Illustration  of  a  typical  spatial  variation  Is  shown  In  Fig.  6. 
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Figure  6.  Spatial  Variation  In  PMP  Intensities 


2.2,2  Temporal  Distribution 

The  factors  governing  the  temporal  distribution  of  the  PHS  are  as 
follows: 

*  PHP  for  all  durations  (up  to  3  days)  occurs  In  the  same  PHS.  The  PHP 
pattern  Is  developed  so  that  any  duration  of  storms  less  than  72  hours  Is 
contained  In  the  PHS. 

*  The  four  6-hr  periods  with  the  greatest  precipitation  may  occur  any  time 
except  during  the  first  24  hours  of  the  storm. 

*  The  6-hr  Increments  of  precipitation  are  arranged  such  that  the 
Increments  decrease  progressively  to  either  side  of  the  greatest  6-hr 
Increment.  An  example  of  one  such  distribution  Is  shown  In  Fig.  7. 

*  The  6-hr  Increments  may  be  distributed  Into  shorter  Intervals.  Fig.  8 
shows  ratios  of  1-hr  to  6-hr  precipitation  for  the  'A'  Isohyet  of  a 
20.000  ml*  storm.  This  ratio  Is  determined  for  the  storm-center 
location  and  used  to  adjust  ratios  read  from  Table  5  for  each  Isohyet 
within  the  storm-area  size.  Haxlmum  5-,  15-,  and  30-m1n  Intervals  are 
given  only  for  the  maximum  1-hr  Increment  within  the  maximum  6-hr 
Increment . 
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Time  {6-hr  Increments) 

1*t  24-hr— 72  hr 
PERIOD 


Figure  7.  Schematic  of  One  Temporal  Sequence  Allowed 
for  6-hr  Increments  of  PMP  (NWS,  198?) 


Figure  8.  Ratio  of  1-hr  to  6-hr  Precipitation  for  'A 
Isohyet  of  a  ?0,000  mi*  Storms  (NWS,  198?) 
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Spatial  Variation  in  PMP  for  4th  through  12th  6-hour  Incromont 
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Set  lion  '< 


PROBABl  I  MAXIMUM  I'Hmi’l  I  Al  ION  DU  I  KMI NA I  I  ON 


3.1  All  Soason  PMP  fsliiiMlps 

Gonerally,  a  throe  slop  process  is  followed  in  dolermininy  I'MI’  in 
rronoroqraphii  regions:  moislurr*  max  Imi/at  ion,  1  ranspos  i  I  ion  and  onvo  lopmonl 
Ihoso  protessos  are  briefly  dost ri bod  below;  HMk  No.  SI  gives  a  (OMiplolc 
explanation  of  the  PHI*  proieduro 

Moisture  maximi/ation  consists  of  increasing  storm  pr ei ipi lat  ion  measured 
in  a  major  hislorli  event  by  a  far  lor  that  reflects  the  maximum  amount  of 
molslun*  that  rould  have  exislei!  in  the  atmosphere  fur  the  storm  Inral  mn  and 
lime  of  ) he  year. 

1 1  ansp'is  i  t  ion  refer-  to  ihc  pfuross  of  moving  a  storm  (i.e.,  isohyol.ii 
palletn)  1 1  oni  Ihr-  lut.ition  wher>-  it  oirtitrr-d  to  anolhor  loialion  ot  inti-resl. 

1  ranspos  i  I  Ion  Is  tart  led  oui  only  wUhin  a  tegion  that  is  homrirp-neous  witti 
rosper  I  to  terrain  and  ini' !  t>or  o  I  tit|y  . 

!  n  VI' lojimen  t  involvi's  loosl.-ui  t  ion  of  -.niorjlh  >  ui  ves  that  i-nve'npe 
pr  ei  ip; '  al  ion  maxima  (or  vat'ieus  I'ui.itions  and  ari-i  ’.i/r”.  'o  oinpe-r.a  I  e  :  oi 
dr.ia  ijafis .  Al-.o  ijeogia()hu  Miimiltiing  is  petlotiin-d  to  ii.'.aM-  ii'ijiinu-i 
t  "OS  I  s  ’..ent  y  . 

Using  Ihosi'  printiples,  PMP's  for  various  sired  area  and  slrtrm  durations 
were  calculated  by  the  NWS.  Ihe  results  are  generalized  and  plotted  in  a 
series  of  figures  (18  4/)  in  HMR  No.  SI.  One  of  those  plots,  the  all  season 
PMP  for  the  6  hr  10-mM  PMP,  is  illustrated  here  in  I  ig.  9. 
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l  itjiiro  All  Season  PMP  (Iruhes)  lor  h  bouts,  ID  mi''  (NWS,  i9/8) 


3.?  Example  Calculation  of  PHP 

It  Is  desired  to  obtain  PMP  estimates  for  the  Leon  River  basin  above 
Belton  Reservoir,  Texas,  fig.  10.  The  drainage  area  Is  3,660  mW;  the 
location  of  the  centroid  of  the  basin  Is  31*45'  N,  98*15*  W. 

From  the  PMP  maps  In  HMR  No.  51.  the  PMP  values  for  area  sizes  larger  and 
smaller  than  the  drainage  area  are  determined.  The  all -season  PMP  for  6  hours 
and  10  m12  was  shown  In  fig.  9  (which  Is  Fig.  18  In  HMR  No.  51).  Values 
taken  from  figures  18  through  47  of  HMR  No.  51  are  shown  In  Table  6. 


Table  6 


PMP  Depths  (by  Area  and  Duration)  for 
Leon  River  Basin 


Area 

(ml2) 

6 

12 

Duration  (hours) 

24 

48 

72 

10 

29.8 

36.2 

41.8 

46.7 

49.8 

200 

22.3 

27.4 

33.0 

37.5 

41.4 

M. 

16.2 

21.2 

26.8 

31.0 

34.5 

9.3 

13.1 

18.1 

22.6 

25.9 

7.2 

10.4 

14.9 

18.8 

21.0 

5.2 

8.2 

11.7 

15.4 

18.4 

From  the  depth-area-duration  data  In  Table  6,  plot  PMP  depth  versus  the 
logarithm  of  drainage  area  for  each  duration.  Draw  smooth  curves  for  each 
duration.  The  curves  should  be  parallel  or  converge  slightly  with  Increasing 
size,  as  Illustrated  In  Fig.  11. 


Figure  11.  Depth- Area-Duration  Curves  for 
Leon  River  Basin  (NWS.  1962) 
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It  Is  highly  recommended  to  physically  plot  the  Initial 
depth-area-duration  values  taken  from  the  various  figures  (figures  18-47)  In 
HHR  51  as  shown  In  Fig.  11  of  this  users  manual.  This  Initial  plotting  of 
the  basic  Input  data  serves  two  functions.  One,  It  eliminates  reader  errors 
from  basic  misinterpretation  of  values  read  off  the  figures  In  HMR  51. 

Second,  Initial  Important  smoothing  of  the  basic  precipitation  data  Is 
applied. 

Using  the  depth-area-duration  graph  of  Fig.  11,  determine  PHP  depths  for 
the  storm-area  sUes  of  Interest.  A  plot  of  the  area-specific  PMP  values 
versus  duration  on  linear  graph  paper  can  be  used  to  obtain  PMP  depth  for  any 
duration  between  6  and  72  hours.  The  above  curve  fitting  Is  done 
automatically  by  the  HMR5?  program.  The  user  need  only  provide  the 
depth-area-duration  data,  e.g.  Table  6,  for  the  geographic  location  1n 
question.  Using  this  PMP  data,  the  PMS  Is  calculated  as  described  In  the 
next  section. 


1/ 


Section  4 


EXAHPLF  CALCULATION  Of  A  PMS 


lo  Illustrate  the  procedures  given  In  HHR  No.  52,  suppose  that  It  Is 
desired  to  derive  a  storm  for  a  storm-area  size  of  3,000  ml^  for  the  Leon 
River  basin.  The  storm  Is  to  be  centered  at  31“45'  N,  98*15'  W  with  an 
orientation  of  314®,  as  shown  In  fig.  12. 


figure  12.  Isohyetal  Pattern  Placed  on  Leon  River  Basin 
(NWS,  1982) 


DEPTH  -  inch** 


4 . 1  Determination  of  PHP  Intensity  versus  Duration 

The  PHP  for  an  area  of  3,000  ml*  (determined  from  Fig.  11)  Is  shown 
In  Fig.  13.  The  Incremental  precipitation  amounts,  calculated  (Including 
additional  smoothing)  from  the  cumulative  amounts  given  In  Fig.  13,  are  as 
follows: 


6-hr 

Increment 


PHP 

Inches 


Increment 


Inches 


IS  24 


42  4B  S4  60  66 


DURATION -hourt 


Figure  13.  PHP  Deplh-Ourallon  for  3,000  ml*  Storm 
on  I  eon  River  Basin 
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4 . 2  Adjustment  for  Storm  Orientation 

From  FJg.  3,  the  preferred  orientation  for  this  storm  location  Is  208®, 
which  differs  by  106®  from  the  selected  orientation  of  314®.  from  fig.  4, 
the  PMP  must  therefore  be  multiplied  by  0.85.  The  adjusted  PHP  Is  as  follows 


6-hr 

PMP 

Increment 

Inches 

Increment 

Inches 

1 

9.69 

7 

0.85 

2 

3.52 

8 

0.85 

3 

2.30 

9 

0.77 

4 

1.79 

10 

0.77 

5 

1.39 

11 

0.68 

6 

1.02 

12 

0.68 

4 . 3  Determination  of  Basin-Average  Precipitation 

The  nomograms  In  Figs.  14  through  17  are  required  for  determining  the 
spatial  distribution  of  precipitation.  Tables  7  through  10  show  computatio 
to  determine  basin-average  precipitation  for  the  four  largest  6-hr 
Increments.  Basin-average  precipitation  for  the  remaining  Increments  (5th 
through  12th}  can  be  obtained  as  a  proportion  of  the  4th  largest  Increment, 
because  the  same  spatial  distribution  Is  used.  For  example,  the  3,000 
ml^  PMP  for  the  5th  largest  Increment  Is  1.39  Inches.  Basin-average 
precipitation  for  the  Increment  Is  1.39  (1.69/1.79)  =  1.31  Inches. 
Basin-average  precipitation  for  the  twelve  Increments  Is  as  follows; 

6-hr  PMP 


Increment 

Inches 

Increment 

Inches 

1 

8.31 

7 

0.80 

2 

3.25 

8 

0.80 

3 

2.17 

9 

0.73 

4 

1 .69 

10 

0.73 

5 

1.31 

11 

0.64 

6 

0.96 

12 

0.64 

4.4  1 emporal  Ar ranqement 

An  acceptable  temporal  arrangement  (Fig.  7)  of  the  basin  average 
precipitation  Is  as  follows: 


6-hr 

Basin  -Average 

6-hr 

Basin-Average 

Da^ 

Per  IwJ 

Precipitation 

Day 

Period 

Precipitation 

1 

1 

0.64  Inches 

7 

1.31 

2 

0.64 

8 

2.17 

3 

0.73 

3 

9 

8.31 

4 

0.73 

10 

3.25 

? 

5 

0.80 

1 1 

1.69 

6 

0.96 

1? 

0.80 
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Table  1 


Basin-Average  Precipitation  for  Largest  6-hr  Increment 


Isohyet 

Area 

Within 

Isohyet 

Area  Within 
Isohyet  and 
Basin  (A) 

Incrmiental 
Area  aA 

X  of  3000 
sq.mi.  PMP 

Isohyet 
Precipi tation 

Average 
Precipitation 
on  aA 

Volume  of 
Precipitation 

sq.mi . 

sq.mi. 

sq.mi . 

(9.69  in.) 

in. 

in. 

in.  -sq.ini . 

A 

10 

10 

10 

191 

18.51 

18.51 

185.1 

B 

25 

25 

15 

179 

17.35 

17.93 

268.9 

C 

50 

50 

25 

166 

16.09 

16.72 

418.0 

0 

100 

100 

50 

154 

14.92 

15.51 

775.5 

E 

175 

175 

75 

142 

13.76 

14.34 

1075.5 

F 

300 

300 

125 

132 

12.79 

13.28 

1660.0 

6 

450 

450 

150 

122 

11.82 

12.31 

1846.5 

H 

700 

700 

250 

112 

10.85 

11.34 

2835.0 

I 

1000 

971 

271 

102 

9.88 

10.37 

2810.0 

J 

1500 

1364 

393 

92 

8.91 

9.39 

3690.3 

K 

2150 

1852 

488 

83 

8.04 

8.48 

4138.2 

L 

3000 

2434 

532 

74 

7.17 

7.61 

4429.0 

M 

4500 

3171 

737 

44 

4.26 

6.01 

4429.4 

N 

6500 

3660 

489 

25 

2.42 

3.80 

1858.2 

30419.9 

Basin  Avg.  Ppt.  =  30419.9  ‘  3660  =  8.31  in. 
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lable  8 

Basin-Average  Precipitation  for  Second  Largest  6 -hr  Increment 


Table  9 


Basin-Average  Precipitation  for  Third  largest  6-hr  Increment 


Isohyct 

Arpa 

Wi thin 
Isohyct 

Ai'o.i  Within 
Isohyct  and 
Basin  (A) 

Incremental 
Area  aA 

— 

%  of  3000 
sq.mi.  PMP 

Isohyet 

Precipitation 

Average 
Precipitation 
on  aA 

Volume  of 
Precipi tation 

sq.mi . 

sq.iiil . 

sq.mi . 

(2.30  in.) 

in. 

in. 

in. -sq.mi . 

A 

10 

10 

10 

105.7 

2.43 

2.43 

24.3 

B 

25 

25 

15 

104.6 

2.41 

2.42 

36.3 

C 

50 

50 

25 

103.5 

2.38 

2.40 

60.0 

0 

100 

100 

50 

102.5 

2.36 

2.37 

118.5 

E 

175 

175 

75 

101.7 

2.34 

2.35 

176.3 

F 

300 

300 

125 

101.3 

2.33 

2.335 

291.9 

G 

450 

450 

150 

100.9 

2.32 

2.325 

348.8 

H 

700 

700 

250 

100.5 

2.31 

2.315 

578.8 

1 

1000 

971 

271 

100.2 

2.30 

2.305 

624.7 

J 

1500 

1364 

393 

99.9 

2.30 

2.30 

903.9 

- 1 

K 

2150 

1052 

408 

99.6 

2.29 

2.295 

1120.0 

L 

3000 

2434 

582 

99.3 

2.28 

2.285 

1329.9 

M 

4500 

3171 

737 

76 

1.75 

2.07 

1525.6 

N 

6500 

3660 

489 

57 

1.31 

1.64 

802.0 

7941.0 

Basin -Avg.  Ppt.  =  7941.0  :  3660  =  2.17  in. 
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lable  10 


Basin -Average  Precipitation  for  Fourth  Largest  6-hr  Increment 


Isohyet 

Area 

Wi  thin 
Isohyet 

Area  Within 
Isohyet  and 
Basin  (A) 

Incremental 
Area  aA 

%  Of  3000 
sq.mi.  PMP 

Isohyet 
Precipi  Lation 

Average 
Precipitation 
on  aA 

Volume  of 
Precipitation 

sq.mi . 

sq.mi . 

sq.mi . 

(1.79  in.) 

in. 

in. 

in.-sq.mi . 

A 

10 

10 

10 

100 

1.79 

1.79 

B 

25 

25 

15 

100 

1.79 

1.79 

50 

50 

25 

100 

1.79 

1.79 

0 

100 

100 

50 

100 

1.79 

1.79 

E 

175 

175 

75 

100 

1.79 

1.79 

limn 

300 

300 

125 

100 

1.79 

1.79 

G 

450 

450 

150 

100 

1.79 

1.79 

H 

700 

700 

250 

100 

1.79 

1.79 

I 

1000 

971 

271 

100 

1.79 

1.79 

J 

1500 

1364 

393 

100 

1.79 

1.79 

K 

2150 

1852 

480 

100 

1.79 

1.79 

L 

3000 

2434 

502 

100 

1.79 

1.79 

M 

4500 

3171 

7.17 

1.36 

1.62 

1193.9 

N 

6500 

3660 

489 

57 

1.02 

_ 

1.28 

625.9 

Basin  Avg.  Ppt.  =  [1.79  (2A3A)  +  1193.9  +  625.9]  :  3660  =  1.69  in. 


AREA 


PERCENT  OP  3<^  6-hr  PMP  INCREMENT  PERCENT  OP  3nd  4-l.r 

PMP  INCREMENT 


Mgure  16.  Nomogram  for  3r(j  6-hr  PMP  Increment 
(NWS,  198?) 
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4 . 5  PMf  De  t e r m  Vna  t » on 

Ihe  example  just  shown  provides  a  storm  for  a  specific  storm-area  s1?e, 
storm  centering  and  storm  orientation.  What  would  be  required  for  a  PMF 
analysis,  however,  would  be  the  storm  that  produces  the  maximum  peak 
discharge  or  runoff  volume,  depending  on  the  project  purpose,  for  the  3,660 
ml*  drainage  basin.  It  Is  therefore  necessary  to  try  various 
combinations  of  storm-area  size,  centering,  orientation  and  temporal 
distribution  until  maximization  of  peak  discharge  or  runoff  volume  Is 
achieved.  Fig.  18  Illustrates  how  basin  average  precipitation  varies  with 
storm-area  size  for  the  Leon  River  basin.  For  the  particular  centering  and 
orientation  used  In  the  example,  a  storm-area  size  of  2,150  ml*  produces 
the  maximum  basin  average  precipitation.  As  may  be  deduced  from  the  example 
problem,  voluminous  computations  are  required  to  determine  the  PMS.  The 
HMR5?  computer  program  has  been  developed  to  facilitate  that  task. 
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figure  18.  18-hr  Has  In -Average  Precipitation  for 

3,660  ml*  Icon  River  Basin 


Section  5 


COMPUUR  PROGRAM  HMR5?  PROCfDURES 


5 . 1  Digital  DefInUlon  of  Basin  Geometry 

The  HMR5?  computer  program  uses  a  digital  definition  of  the  watershed 
boundaries  for  computing  basin-average  precipitation  from  watershed  area  and 
superposed  Isohyetal  patterns.  Ihe  boundary  of  a  drainage  basin  Is  defined  by 
line  segments  Joining  a  sequence  of  coordinate  points.  The  sequence  of 
boundary  points  should  be  counter-clockwise  around  the  basin  or  subbasin.  If 
the  direction  Is  clockwise,  the  program  will  reverse  the  order  to  be 
counter-clockwise  for  later  calculations. 

5.1.1  Georoe  t  r  1  c_Proper  1 1  e  s 

The  following  geometric  properties  of  the  watershed  are  calculated  using 
the  n  boundary  points  having  coordinates  (x^,  y^): 

Area  (A): 


A  ^  -  (1/?)  I  (x^^l  -  x^)  f  y^) 

U1 


(1) 


Centroid  coordinates  [%(■,  y^)  ; 


Xc  =  -  [1/(6A)]  I  ^  -  x,)  f  y,^^  y,  f  y\) 


.  .  (?) 


Yc  ^  '  ni/hA)]  X  ^  -  x^)  (y^^^  *'  ’‘i  > 

^  Y)  .... 


■  .  (3) 


Moment  of  Inertia  about  x  and  y  axes  (I^,  ly)  : 

n 


1  ^  -  (1/1?)  i  (X, 

J  \ 


t  ?  V  <  ' 

1  + 1  ■  ■  ■ 


Itl 

«,)  (y^.,  •  y^)  (y,.,  *  y|)  •  . 

•  .  (4) 

i^l 

■•1^  •  ^  ^  Y,) 

w 

■'  t  H 

y,)  ^  x/  (y,^,  .  3y,)]  .  .  .  . 

•  .  (5) 
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Product,  of_  inertia  about  the  orlgtn  (P  )*• 

*  “•  X  y 

Pxy  =  -  (1/24)  I  ^  -  xp  ^  *1^ 

"  y^  <*U1  "  . 

Angle  to  rotate  coordinate  axes  to  produce  minimum  moment  of  Inertia  about 
the  x-axIs. 

(e^)  =  (1/2)  arctan  [ -^P^y/^x  ^ 

(6^)  .  (1/2)  arctan  {[ -PP^^/d^-Iy)]  W21  for  Ix  >  ly  •  •  • 

5.1.2  C^r dinate  Systems 

Ihe  basin  Is  described  In  an  x,y  coordinate  system  with  x-axIs  directed 
eastward  and  y-axis  directed  northward.  Placement  of  the  coordinate  system 
origin  and  the  coordinate  units  are  arbitrary. 

The  Isohyetal  pattern  of  the  PHS  Is  described  In  a  u,v  coordinate  system 
which  has  Its  origin  at  the  storm  center  and  axes  parallel  to  the  major  and 
minor  axes  of  the  elliptical  pattern.  Ihe  coordinate  units  are  In  miles. 

The  transformation  from  the  x.y  axes  to  the  u,v  axes  Is  given  by: 


u  =  (X  -  Xj)  cos  e  t  (y  -  ys)  sin  e  Fj . (8a) 

V  =  (y  -  Ys)  e  -  (x  -  Xj)  sin  0  F5 . (0b) 


where  (Xj.yj)  Is  the  position  of  the  storm  center,  0  Is  the  angle  of 
rotation  from  the  x-axIs  to  the  u-axis,  and  F^  Is  a  scale  factor  In  miles 
per  x.y  coordinate  unit. 

5 . 2  Cj^lcujadon  of^  PT1P  f or  a_Storm  Area 

This  calculation  Is  made  for  a  given  storm  area  and  depth -area -duration 
data  from  HMR  No.  51.  The  procedure  used  Is  as  follows: 

5.2.1  Interpolation  of  PAP_()eptji-Area  Duration  Cury_es 

The  user  supplies  PMP  depths  from  HMR  No.  51  for  standard  areas  and 
durations  that  bracket  the  storm  area  of  Interest.  Straight  line 
Interpolation  Is  used  In  the  program  to  calculate  PMI’  for  the  storm  area  (or 
durations  of  6,  12,  24,  48,  and  7?  hours  (e.g.  fig.  19). 
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Between  10  and  200  ml^.  the  error  from  using  straight  line 
Interpolation  versus  curve  fitting  Is  significant.  To  reduce  this  error, 
additional  point  Is  calculated  midway  between  log]Q  10  and  log^g  200 
ro1».  This  area  Is  approximately  44.7  m1».  Fig.  19.  The  difference  In 
PMP  from  10  to  200  ml*  for  this  curve  was  8.5  Inches.  To  estimate  PHP  at 
44.7  ml*  for  other  locations  an  adjustment  of  1/8.5  times  the  difference 
In  PMP  from  10  to  200  ml*  Is  added  to  the  PMP  resulting  from  straight  line 
Interpolation. 


DUMATMM 


5.2.2  PMP  Depth-Duration  Curve  Fitting 

A  logarithmic  curve  Is  fit  through  the  depth  versus  log  duration  data 
using  a  least-squared-error  fit. 

The  equation  of  the  curve  Is 

P  =  a  t  b  in(DtS) .  (9) 

where  P  Is  PHP  In  Inches;  0  Is  duration  In  6-hr  periods;  S  Is  a  shift 
adjustment  ranging  from  -1  to  2.5  6-hr  periods;  and  a  and  b  are  derived 
constants.  The  shift  adjustment  Is  determined  by  trial  and  error  to  minimize 
the  sum  of  squared  errors  to  reproduce  a  curve  such  as  was  shown  In  Fig.  13. 

5.2.3  Incremental  PHP  Calculation 

PMP  Is  computed  from  equation  (9)  for  durations  from  6  to  72  hours. 
Incremental  6-hr  PMP  Is  then  computed  as  the  difference  between  the  values  of 
PMP  for  successive  durations.  The  PMP  for  the  1*^  6-hr  period,  P^,  Is 

P^  =  a  +  b  in(l  >  s)  -  P^-1 .  (10) 


5 . 3  Calculation  of  PMS  for  Given  Storm  Area  and  Orientation 
5.3.1  Orientation  Adjustment 

PMP  Is  calculated  for  the  storm-area  size  as  described  In  Section  3. 

This  PMP  Is  multiplied  by  the  orientation  factor,  f^,  which  Is  computed  from 

fo  =  1  -  0.15  Cl  Cp .  (11) 


where  Ci^  < 


Oc  - 


Op  -  40 


25 

1 


.  -  ,  A.  -  300 

and  Cp  =  <  _ 

'  2700 

1 


Is  the  storm  area. 


for 

|0s  - 

OpI 

<  40 

for 

40  < 

10s 

-  OpI  <  65 

for 

lOs  - 

OpI 

>  65 

for 

^s  5 

300 

for 

300  < 

^s 

<  3000 

for 

*S  - 

3000 

Ion, 

Op  Is 

the 

preferred 
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(12) 


Ihe  adjusted  PMP,  P^',  for  the  6-hr  period  Is  calculated  as 

Pr  ^  Pi  *  Fo . 


5.3.2  Isohyet  values 

Ihe  percent  of  PMP  for  each  Isohyet  Is  Interpolated  from  the  data  In 
Tables  1  through  4  using  straight  line  Interpolation  of  percent  versus  natural 
log  of  area.  Hultiplying  PMP  by  the  percentages  gives  the  precipitation 
amount  for  each  Isohyet  for  each  6-hr  period  for  a  given  storm-area  size. 

Thus,  the  PMP  for  the  J^F*  Isohyet  and  I^Fi  6-hr  period  Is 

Pl.j  ^  pr  ‘  pcij . (13) 

where  PCTj  Is  the  percentage  of  the  PMP  In  the  Isohyet. 

5.3.3  Calculation  of  area  enclosed  bjL aiL ellipse  and  t h^ basin  boundary 

Coordinates  of  the  boundary  points  are  transformed  to  the  PMS  coordinate 
system  using  equations  (8a)  and  {8b).  The  Intersections  of  the  ellipse  with 
the  basin  boundary  are  located.  Each  arc  of  the  ellipse  between  Intersections 
Is  approximated  by  20  points  located  at  equiangular  Increments.  The  area  of 
the  basin  encompassed  by  the  ellipse  (see  Flq.  20)  Is  calculated  from  equation 
(1). 


w 


figure  ?0,  basin  Area  Encompassed  by  Ellipse 
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5.3.4  B^Jn  a^eracie  precipitation 

The  precipitation  volume,  V,  between  Isohyets  Is  computed  as  the  volume 
of  a  truncated  cone  (see  Tig.  ?1).  The  volume  for  the  l*^”  6-hr  period  Is 

V,  .  (Aj  .  AjAj,,  .Aj.,) . 

where  subscript  J  Identifies  the  Isohyet  and  A  Is  the  drainage  basin  area 
encompassed  by  the  ellipse  corresponding  to  the  Isohyet.  Basin  average 
precipitation  Is  the  sum  of  volumes  between  Isohyets  divided  by  the  drainage 
basin  area. 


figure  ?1.  Representation  of  Precipitation  Volume  Between 
Iwo  Isohyets  with  a  Irunrated  Cone 


5.4  Calculation  of  PMS  for  L'me  Jn tergal s  Less  Than  6_  Hours 

After  the  storm-area  size  and  orientation  have  been  determined,  the 
program  calculates  the  temporal  distribution.  The  required  data  for  this 
calculation  are:  the  desired  time  Interval,  at.  In  minutes;  and  the  ratio 
of  1 -hr  to  6-hr  precipitation  for  the  20,000  ml*  'A'  Isohyel,  referred  to 
as  R16A20,  from  Fig  8.  The  procedure  used  In  the  HMR52  program  Is  as  follows. 

The  NWS  has  developed  a  set  of  curves  describing  the  1-hr/6-hr  ratio  for 
each  Isohyet  and  various  storm-area  sizes  (figure  40,  NWS,  1982).  Those 
curves  have  been  tabulated  by  HEC  Into  a  set  of  discrete  values  as  shown  In 
Table  5.  Those  Table  5  values  are  automatically  specified  In  the  HMR5? 
program  (see  Computer  Requirements  section  and  Table  24).  The  value  of  that 
ratio  for  the  Isohyel,  RlblABj,  for  a  specific  storm-area  size.  Is 
Interpolated  from  Table  5  data  using  straight  line  Interpolation  of  logarithms 
of  the  area.  Those  generalized  values  are  adjusted  bv  the  required  Input 
value,  RlbA20,  to  obtain  the  required  ratio  for  the  J*^  Isohyel,  Rlbj- 

R16j  .  RlbTABj  t  R16A20 . (15) 

Ihe  1 -hr  rainfall  for  the  J^*’  Isohyet  Is  computed  as 

^’1J  =  Rl&J  *  Pl,J . (16) 

where  Pi  i  Is  the  maximum  6-hr  precpUatlon  for  Isohyet  J  calculated  In 
equation’ (13) . 

A  least  square  quadratic  equation  Is  fit  through  the  0-,  1-,  6-,  12-, 

18-,  and  24-hr  precipitation  amounts  for  each  Isohyet.  For  residual 
precipitation  (outside  the  storm  area)  the  1-hr  amount  Is  omitted. 

Precipitation  Is  Interpolated  for  5-,  10-,  15-,  30-m1n,  and  3-hr  durations 
for  each  Isohyet.  The  average  precipitation  over  the  basin  Is  computed  for 
each  of  these  durations.  Fig.  22. 


Owr«f(«ii 


I  pjre  22.  PMS  Depth  versus  Duration 
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Incremental  precipitation  Is  calculated  for  the  total  number  of  At 
Intervals  from  the  basin -averaged  depth  versus  duration  data.  This  results  In 
a  sequence  of  Intervals  with  decreasing  precipitation  Intensity,  Fig.  23. 
Incremental  precipitation  Is  assumed  to  be  uniform  within  each  6-hr  period 
beyond  ?4-hrs  duration.  These  Intervals  are  rearranged  to  form  the  PMS  as 
follows : 


The  position  of  the  largest  6-hr  Increment  may  occur  any  time  after 
the  first  ?4  hours  of  the  storm  as  shown  In  fig.  24,  The  seventh 
position  (hours  37  42)  Is  chosen  by  default. 

Ihe  largest  At  Increment  of  precipitation  Is  placed  In  the  middle 
of  the  largest  6hr  precipitation  Interval.  The  remaining  At 
Increments  for  that  6-hr  period  are  arranged  alternately  before  and 
after  the  largest  Increment. 

The  remaining  6-hr  Intervals,  with  decreasing  precipitation 
magnitude,  are  arranged  alternately  before  and  after  the  largest 
6-hr-prec1p1 tatlon  Interval,  except  the  second,  third,  and  fourth 
largest  6-hr  Increments  cannot  be  placed  In  the  first  24  hours  of  the 
storm.  Ihese  Increments  are  placed  after  the  largest  Increment,  If 
their  normal  position  would  fall  In  the  first  24  hours  of  the  PMS. 

The  At  Increments  within  the  second,  third  and  fourth  largest  6-hr 
Increment  are  arranged  to  Increase  towards  the  storm's  peak 
precipitation.  This  results  In  a  trlangular-1 ike  hyetograph  of  PMS 
precipitation. 
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figure  23,  Incremental  PMS  Depth  Histogram 


Section  6 


OPTIMIZATION  Of  STORM-APEA  SIZE  AND  ORIENTATION 


The  amount  of  precipitation  on  a  basin  Is  affected  by  the  storm 
placement,  storm-area  size  and  storm  orientation.  HMRS?  uses  a  piocedure  to 
estimate  storm-area  size  and  orientation  which  will  produce  maximum 
precipitation  on  the  basin.  This  procedure  will  determine  the  optimal 
storm-area  size  and  orientation  for  most  basins.  However,  because  of  the 
interaction  of  storm  placement  and  orientation,  several  trials  should  be  made 
to  verify  that  the  optimal  values  have  been  found. 

6 . 1  Storm  Center 

When  the  storm  placement  Is  not  given,  the  storm  center  Is  placed  at  the 
basin  centroid. 

6 • 2  Storm-Area  Size 

If  the  storm  orientation  Is  not  given,  a  trial  orientation  along  the  axis 
for  which  the  basin  has  a  minimum  moment  of  Inertia  Is  used.  This  orientation 
Is  selected  as  an  analytically  determinable  orientation  which  Is  most  likely 
to  produce  maximum  precipitation  on  the  basin  without  regard  to  the 
orientation  adjustment  factor.  Using  the  trial  orientation  the  average 
precipitation  on  the  basin  Is  calculated  for  several  storm-area  sizes  (see 
example  1s  Section  9).  The  storm-area  size  which  produces  maximum  precipita¬ 
tion  Is  selected  as  the  critical  storm-area  size.  The  critical  storm-area 
size  may  need  to  be  recomputed  If  other  orientations  or  placements  are  used. 

6 . 3  Storm  Orientation 

Because  of  Interaction  between  basin  shape  and  the  orientation  adjustment 
factor,  the  trial  orientation  used  to  select  storm-area  size  may  not  produce 
maximum  precipitation  on  the  basin.  Using  the  critical  size  selected  above, 
the  storm  precipitation  Is  calculated  for  orientations  at  10-degree  Increments 
between  135  degrees  and  315  degrees.  Storm  precipitation  Is  then  calculated 
for  orientations  at  plus  or  minus  5  degrees  from  the  previous  best 
orientation.  The  orientation  which  yields  maximum  precipitation  on  the  basin 
Is  chosen  as  the  critical  orientation. 

If  the  critical  orientation  Is  not  the  same  as  was  used  to  select  the 
critical  storm-area  size,  the  storm-area  size  may  have  to  be  recomputed  using 
the  critical  orientation.  This  Is  the  user's  responsibility.  It  Is  not  done 
automatically. 

6 . 4  User  Control  of  Optimization 

Precipitation  on  the  basin  for  the  purposes  of  selecting  storm-area  size 
and  orientation  Is  the  cumulative  precipitation  for  the  3  largest  6-hr 
periods.  Ihe  number  of  periods  may  be  changed  by  the  user. 

Storm  placement,  storm-area  size,  or  orientation  may  be  fixed  by  the 
user.  In  addition  the  trial  orientation  used  to  select  critical  storm-area 
size  may  be  specified  by  the  user. 
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Section  7 


INPUT  DATA  REQUIREMENTS 


Input  data  are  described  In  detail  In  the  Appendix.  Table  11  shows  the 
sections  of  an  Input  data  file  required  to  calculate  a  PMS. 

Oralnage-Basin-Georoetry  data  Includes  coordinates  of  points  on  the  basin 
boundary  and  a  scale  factor. 

Required  Hydrometeorological  data  are  the  preferred  storm  orientation 
from  F1g.  3  and  PMP  estimates  from  figures  18-47  In  HHR  No.  51. 

The  Storm-Specification  data  define  the  storm-area  size,  orientation,  and 
location  of  the  storm  center.  Calculation  of  a  temporal  storm  distribution 
requires  the  desired  time  Interval  and  ratio  of  1-hr  to  6-hr  precipitation 
from  Fig.  8. 

Three  example  Input  data  sets  are  given  In  this  document: 

Leon  R1ver--1ab1e  11 
Jones  Reservo1r--Tab1e  14 
Ouachita  R1ver--Append1x  A  Example 

The  Leon  River  application  Is  for  computation  of  the  PMS  over  a  single 
basin  for  given  storm-area  size,  orientation,  centering,  and  time  pattern. 

The  Jones  Reservoir  application  Is  for  a  multl-subbasin  river  basin  In 
which  the  storm-area  size  and  orientation  are  optimized  by  the  program  for  the 
river  basin  as  a  whole.  The  storm  centering  was  not  specified  so  the  default 
centering  of  the  basin  centroid  Is  used.  The  time  pattern  Is  given.  After 
computation  of  the  PMS  for  the  entire  rive  basin,  the  PMS  1s  calculated  for 
each  subbasin. 

The  Ouachita  River  example  Is  similar  to  the  Jones  Reservoir  application 
except  that  given  storm-area  size,  orientation  and  centering  data  are  provided 
after  the  optimization  Is  completed.  That  given  PMS  Is  then  calculated  and 
subbasin  precipitation  Is  computed  for  the  given  PMS  data,  not  the  optimized 
values. 
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Table  11 


Sample  Input  for  HNR5? 


Job  Oescrlotlon  (Uentifloatl 

TO  PROBA^r,?:  MAXIMUM  STORM  CALCULATION  FOR  KXAMPLR  \  TN  MMR  NO. 
ID  LEON  RI’TER  AT"  BELTON  RSSER'/OtR 


Oralfiaqo  Basil  Ceomet'si 


BN 

LEON 

BS1.0062 

8X 

4^5.4 

53.4 

52.2 

52.5 

46.1 

42.8 

34.3 

30.0 

BX 

13.7 

16.5 

14.2 

12.5 

2.1 

-1.2 

-1.5 

-3.9 

BX 

-12.4 

-19,3 

-20.7 

-25.3 

-28.1 

-30.5 

-31.3 

-38.5 

BX 

-Bl.*) 

-53,5 

-55.1 

-51.5 

-45.5 

-36.9 

-36.7 

-29.5 

BX 

-25. B 

-22.5 

-13.5 

-9.7 

-4.7 

-2.5 

3.5 

5.5 

BX 

5.B 

9.5 

13.1 

15.5 

19.8 

23.7 

26.5 

29.4 

BX 

44.0 

BY 

-45.9 

-33.4 

-3-2.3 

-29.7 

-27.7 

-19.7 

-15.5 

-10.7 

BY 

-1.2 

3.4 

5.5 

3.5 

9.5 

15.9 

19.5 

22.3 

BY 

42.0 

41. B 

44.9 

48.7 

55.4 

56.4 

50.1 

49.3 

BY 

42.1 

33.5 

31.4 

28.3 

28.2 

14.5 

10.0 

10.5 

BY 

-1.5 

-2,9 

-9.5 

-7.8 

-9.4 

-14.0 

-14.5 

-17.2 

BY 

-23 . 5 

-34.0 

-15.7 

-41.0 

-43.4 

-46.7 

-45.5 

-41.5 

BY 

-46.9 

HyHcometiaorolocf  i.7al 

bata  from 

HMR  No. 

51  and 

HO 

203 

HP 

10 

29.3 

36.2 

41.3 

46.7 

49.3 

HP 

200 

22.3 

27.4 

3,3.0 

37.5 

41.4 

HP 

1000 

15.2 

21.2 

25.3 

31.0 

34.5 

HP 

5000 

9.3 

13.1 

18.1 

22.5 

25.9 

HP 

10000 

7.2 

10.4 

14.9 

13.8 

21.0 

HP 

20000 

5.2 

3.2 

11.7 

15.4 

13.4 

Btorm 

Bnec ' 

f  lcatli7n 

BA 

3000 

134 

BC 

-■)  7 

2.4 

BT 

120 

306 

24.7 

22.2 

-7.3 

-7.7 

-42.4 

-44.8 

-29.1 

-25.5 

3.2 

3.6 

35,4 

38.9 

-8.0 

-2.5 

33.5 

37.7 

47.5 

43.9 

7.5 

4.5 

-20.7 

-22.5 

-41.0 

-42.0 

End-of-Joh 


ZZ 


■» 
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Section  8 
PROGRAM  OU1PU1 


8  -1  Pi  Vo. tout 

HMR52  printed  output  begins  with  a  list  of  the  Input  data.  Table  12(a). 

If  data  were  entered  using  the  free-format  option,  this  list  shows  the  data 
values  In  their  proper  fields. 

Table  12(b)  shows  Input  PMP  from  HMR  No.  51  and  PHP  Increments  for  each 
6-hr  Interval  and  <;tandard  storm-area  size  as  Interpolated  from  the  HMR  No.  51 
depth-area -duration  data. 

Table  12(c)  shows  coordinates  of  subbasin  boundary  points  and  subbasin 
area  and  centroid  location  calculated  from  these  coordinates. 

Table  12(d)  shows  the  basin  area  within  each  Isohyet  and  the 
precipitation  amount  assigned  to  the  Isohyels  for  each  6-hr  Interval.  This 
table  also  shows  the  basin-average  depths  calculated  from  these  areas  and 
Isohyet  values. 

When  a  temporal  distribution  Is  requested  for  an  Interval  less  than  6 
hours,  HMR52  computes  a  depth  versus  duration  relation  for  each  Isohyet. 

These  relations  are  shown  In  Table  12(e).  They  are  used  to  calculate  an 
average  depth  versus  duration  relation  which  Is  used  to  calculate  Incremental 
precipitation  for  the  temporal  distribution  In  Table  12(f). 

8.2  Error  Messages 

Ihe  HMR52  program  recognizes  some  Input  and  computational  errors  and 
prints  error  messages  accordingly.  These  error  messages  have  an 
Identification  Number  and  Title;  Appendix  B  contains  an  explanation  of  each 
error  message. 

When  an  error  Is  detected  by  HMR5?,  the  program  will  read  through  the 
remaining  data  and  check  for  Input  errors.  No  calculations  will  be  made  using 
the  remaining  data. 

Ihe  computer  operating  system  may  also  print  error  messages.  When  an 
error  occurs,  the  user  should  first  ascertain  If  It  Is  generated  by  HMR52  or 
by  the  system.  If  It  Is  generated  by  HMR52  (I.e.,  In  the  format  given  In 
Appendix  B)  refer  to  that  appendix  and  take  the  Indicated  actions.  If  the 
error  Is  system  generated,  computer  systems  personnel  should  be  contacted  to 
ascertain  the  meaning  of  the  error.  If  these  system  errors  cannot  be  resolved 
In-house  or  If  there  Is  an  error  In  the  HMR52  program,  the  HEC  should  be 
contacted. 
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Table  1? 


Sample  Output  from  HMR5? 


TAflLR  12(a) 

Sample  Oiitixit  from  HMR52 
HBC  PROAARLR  MAXIMTIM  STORM  INPOT  DATA 


MR 

Tn  . 

_ 1  .  ^ 

_ 5  ^ 

_ - 

. 4.  . 

. A.. 

. 6.  a 

. 7.. 

. 8. . 

. 9. 

....10 

i 

ID 

PROBABLE  MAXIMUM  STORM 

calculation  POR 

EXAMPLE 

1  IN  HMR 

NO.  52 

2 

in 

LEOM  RTTER  AT  BELTON  RESERVOIR 

3 

BN 

LEON 

s 

BX 

46.4 

53.4 

52.2 

52.5 

46.1 

42.8 

34.3 

30.0 

24.7 

22.2 

6 

BX 

18.7 

16.6 

14.2 

12.6 

2.1 

-1.2 

-1.5 

-3.9 

-7.3 

-7.7 

7 

8X 

-12.4 

-19.3 

-20.7 

-25.8 

-28.1 

-10.5 

-31.1 

-38.6 

-42.4 

-44.8 

8 

ex 

-51.5 

-53.6 

-55.1 

-51.6 

-45.6 

-36.9 

-15.7 

-29.5 

-29.1 

-25.6 

9 

BX 

-25.8 

-22.6 

-18.6 

-9.7 

-4.7 

-2.5 

1.6 

5.5 

1.2 

3.6 

10 

ex 

5.8 

9.6 

13.1 

15.5 

19.8 

-M.7 

26.6 

29.4 

35.4 

18.0 

11 

ex 

44.0 

12 

BY 

-45.9 

-38.4 

-32.8 

-29.7 

-27.7 

-19.7 

-15.5 

-10.7 

-8.0 

-2.6 

11 

BY 

-1.2 

3.4 

5.5 

8.5 

9.5 

^  S.9 

19.5 

22.1 

31.5 

17  • 

14 

BY 

42.0 

41.8 

44.9 

48.7 

55.4 

56.4 

50.1 

49.8 

4"'. 5 

41  4 

15 

BY 

42.1 

13.6 

11.4 

28.8 

28.2 

14.6 

10.0 

*»  .  f. 

l« 

BY 

-1.6 

-2.9 

-9.6 

-7.8 

-9.4 

14.0 

'  1  4  .  ' 

-w  ? 

.  • 

17 

BY 

-23.5 

-34.0 

-36.7 

-41.0 

-41.4 

•  4A  .  T 

-  4  .  ‘'i 

-4'  . 

i 

18 

BY 

-46.9 

*« 

19 

HO 

208 

20 

HP 

10 

29.8 

36.2 

41.8 

4<.  ■» 

1'4  ,  Q 

21 

HP 

200 

22.1 

27.4 

.0 

4 

22 

HP 

1000 

V6.2 

21.2 

26.8 

'K « 

21 

HP 

5000 

9.1 

13.1 

l«.l 

24 

HP 

lOOOO 

7.2 

10.4 

25 

HP 

20000 

5.2 

s,  ? 

11.- 

• ' ,  i 

26 

HA 

1000 

114 

27 

sc 

-2.2 

2.4 

28 

ST 

120 

.106 

29 

22 
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!************************************** 


PROBABLE  MAXIMUM  STORM  CALCULATION  FOR  EXAMPLE  1  IN  HHR  NO.  52 

LEON  RIVER  AT  BELTON  RESERVOIR 

PROGRAM  IS  GIVEN  THE  STORM  AREA  AND  ORIENTATION 


PMP  DEPTHS  FROM  HMR  51 


AREA 

(SO.  MI.) 

6-HR 

12-HR 

DURATION 

24 -HR 

48-HR 

72 -HR 

10. 

29.80 

36.20 

41.80 

46.70 

49.80 

200. 

22.30 

27.40 

33.00 

37.50 

41.40 

1000. 

16.20 

21.20 

26.80 

31.00 

34.50 

5000. 

9.30 

13.10 

18.10 

22.60 

25.90 

10000. 

7.20 

10.40 

14.90 

18.80 

21.00 

20000. 

5.20 

8.20 

11.70 

15.40 

18.40 

STORM  AREA 

PMP 

DEPTHS 

FOR  6- 

HOUR  INCREMENTS 

10. 

29.71 

6.47 

3.29 

2.22 

1.68 

1.35 

1.13 

.97 

.85 

.76 

.68 

.62 

25. 

27.96 

6.17 

3.25 

2.22 

1.69 

1.36 

1.14 

.98 

.86 

.77 

.69 

.63 

50. 

26.51 

5.91 

3.21 

2.21 

1.69 

1.37 

1.15 

.99 

.87 

.78 

.70 

.64 

100. 

24.37 

5.54 

3.14 

2.21 

1.70 

1.38 

1.17 

1.01 

.89 

.79 

.72 

.65 

175. 

22.65 

5.25 

3.09 

2.20 

1.71 

1.40 

1.18 

1.02 

.90 

.81 

.73 

.67 

300. 

20.70 

5.18 

3.05 

2.17 

1.69 

1.38 

1.17 

1.01 

.89 

.80 

.72 

.66 

450. 

19.14 

5.21 

3.03 

2.14 

1.66 

1.36 

1.15 

.99 

.87 

.78 

.71 

.65 

700. 

17.46 

5.21 

3.00 

2.12 

1.64 

1.33 

1.13 

.97 

.86 

.77 

.69 

.63 

1000. 

16.10 

5.23 

2.98 

2.09 

1.62 

1.32 

1.11 

.96 

.85 

.76 

.68 

.62 

1500. 

14.37 

4.91 

2.90 

2.07 

1.61 

1.31 

1.11 

.96 

.85 

.76 

.69 

.63 

2150. 

12.83 

4.64 

2.83 

2.04 

1.59 

1.31 

1.11 

.97 

.85 

.77 

.69 

.63 

3000. 

11.40 

4.39 

2.75 

2.01 

1.59 

1.31 

1.11 

.97 

.86 

.77 

.70 

.64 

4500. 

9.67 

4.09 

2.66 

1.98 

1.57 

1.31 

1.12 

.98 

.87 

.78 

.71 

.65 

6500. 

8.39 

3.88 

2.50 

1.85 

1.47 

1.22 

1.04 

.91 

.81 

.72 

.66 

.60 

10000. 

7.06 

3.64 

2.27 

1.65 

1.30 

1.07 

.91 

.80 

.70 

.63 

.57 

.52 

15000. 

5.98 

3.18 

2.12 

1.60 

1.28 

1.07 

.92 

.80 

.71 

.64 

.58 

.54 

20000. 

5.23 

2.84 

2.00 

1.55 

1.26 

1.07 

.92 

.81 

.73 

.66 

.60 

.55 
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A. 


TAn’.R  l’(c) 

Sample  Ontniit  Trom  WMRS?  (contlnup'l) 


BOUNDARY  COORDINATES  FOR  IRON 


X 

4S.4 

S3. 4 

S2.7 

S7.S 

46.1 

42.9 

34.3 

30.0 

24.7 

22.7 

Y 

-4S.9 

-3(3.4 

-32.9 

-29.7 

-27.7 

-19.7 

-IS.S 

-10.7 

-8.0 

-2.6 

X 

IB. 7 

Ifi.B 

14.7 

12.6 

2.1 

-1.7 

-l.S 

-3.9 

-7.3 

-7.7 

Y 

-1.7 

3.4 

S.S 

S.S 

9.S 

IS. 9 

19.  S 

22.3 

33. S 

37.7 

X 

-17.4 

-19.3 

-20.7 

-2S.9 

-28.1 

-30.  S 

-31.3 

-38.6 

-42.4 

-44.8 

Y 

42.0 

41.9 

44.9 

49.7 

SS  4 

S6.4 

SO.l 

49.9 

47. S 

43.9 

X 

-SI  .A 

-S3.S 

-SS.l 

-SI.  6 

-4S.6 

-36.9 

-36.7 

-29. S 

-29.1 

-25.6 

Y 

47.1 

33.6 

31.4 

78.9 

78.2 

14.6 

10.0 

10.6 

7.6 

4.6 

X 

-7S.B 

-27.6 

-19.6 

-9.7 

-4.7 

-2.S 

3.6 

S.S 

3.7 

3.6 

V 

-1  .6 

-2.9 

-9.6 

-7.9 

-9.4 

-14.0 

-14. S 

-17.2 

-20.7 

-22,5 

X 

S.B 

9.6 

13.1 

IS.S 

19.9 

23.7 

26.6 

20.4 

35.4 

38.9 

Y 

-73. S 

-34.0 

-36.7 

-4t  .0 

-43.4 

-46.7 

-4S.6 

-4'  .6 

-41.0 

-42.0 

X  44  .n 
Y  -4<;.9 


SCAIE  =  l.nofil  Mir,F,S  PER  COORDINATE  UNIT 
BASIN  AREA  *  ISfiO.S  SO.  MI. 

BASIN  CENTROID  COORDINATES,  X  =  Y  «  7.4 


TAOLR  12(^1 

Samol«»  Otitnut  from  HMR52  (continne'l) 


ARKA 


T 

3  OR  YET 

WITHIN 

ARPA 

BASIN 

(: 

I'l.MI.I 

ISO. Ml. I 

A 

10. 

10. 

n 

26. 

26. 

c 

60. 

60. 

D 

100. 

100. 

E 

176. 

176. 

F 

300. 

100. 

G 

460. 

460. 

H 

700. 

687. 

I 

1000. 

976. 

J 

1500. 

1397. 

K 

2160. 

1888. 

L 

3000, 

2442. 

M 

4600. 

3246. 

N 

6600. 

3661. 

0 

lOOOO. 

3660. 

P 

15000. 

1660. 

Q 

26000. 

1660. 

n 

40000. 

3660. 

60000. 

3660. 

STORM  AREA 


18. SI 
'7.34 
'S.09 
14.92 
'3.76 
12.79 
11.82 
10. 8S 
9.88 
8.9V 
8.04 
7.17 

4.26 
2,42 

1.26 
0.39 
0.00 
0.00 
0.00 


PRORART.E  MAXIMUM  STORM  FOR  LI^nM 
3000.  SQ.  MI.,  ORIENTATION  =  134.,  PREFERREtl  ORIPNTATTON 
STORM  CENTER  COOROINATES,  X=  -2.2,  v=  2.4 


=  2''3. 


AVERAOE  OEPTH 


8.24 


2 

4.46 

4.33 

4.20 

4.11 

4.03 

3.96 

3.88 

1.81 

3.78 

1.70 

3.62 

1.68 

2.60 

1.68 

1.06 

0.14 

0.00 

0.00 

0.00 

3.41 


2.47 

2.46 

2.42 

2.40 

2.38 

2.37 

2.16 

2.35 

2.36 
2.34 
2.13 
2.32 
1.78 
V.33 
0.88 
0.28 
0,00 
0.00 
0.00 

2.19 


1.71 

1.71 

1.71 

1.71 

1.71 

1.71 

1.71 

1.71 

1.71 

1.71 

1.71 

1.71 

1.30 

0.97 

0.64 

0.21 

0.00 

0.00 

0.00 

1.60 


1.36 

1.35 

1.36 
1.35 

1.35 

1.36 
1.36 

1.35 

1.36 
1.35 
1.35 
1.35 
1.02 
0.77 
0.51 
0.16 
0.00 
0.00 
0.00 

1.26 


l.ll 

1.11 

1.11 
1.11 
l.ll 
l.ll 
l.ll 
l.ll 
l.ll 
l.ll 
l.ll 
l.ll 
0.85 
0.63 
0.42 
0.13 
0.00 
0.00 
0.00 

1.04 


lOR  1NCREM13NTS 
7  8 

OP  PMS 

9 

'  1 

11 

17. 

0.95 

0.82 

0.73 

0.66 

0.59 

0.54 

0.96 

0.82 

0.73 

0.66 

0.69 

0.54 

0.96 

0.82 

0.73 

0.66 

0.59 

0.54 

0.96 

0.82 

0,73 

0,66 

0.59 

0.54 

0.96 

0.82 

0.73 

0,66 

0.59 

0.54 

0.95 

0.82 

0.73 

0.66 

0.59 

0.54 

0.95 

0.82 

0.73 

0.66 

0.59 

0.54 

0.95 

0.82 

0.73 

0.66 

0,59 

0.54 

0.96 

0.8.’ 

0.71 

0.66 

0.59 

0.54 

0.96 

0.82 

0.73 

0.66 

0.59 

0.54 

0.96 

0.82 

0.73 

0.66 

0.59 

0.54 

0.95 

0.82 

0.71 

0.66 

0.59 

0.54 

0.72 

0.63 

0.56 

0.50 

0.45 

0.41 

0.54 

0.47 

0.42 

0.37 

0.34 

0.51 

0.36 

0,31 

0.27 

0,25 

0.22 

0.20 

0.11 

O.IO 

0.09 

0.08 

0.07 

0.07 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0,00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.89 

0.77 

0.68 

0.61 

0.56 

0.51 
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TABLE  l?fAl 

S^ifiplA  Output  from  HMR5?  (contin^od) 


TIME  INTERVAL  «  120.  MINUTES 

1-HR  TO  fi-MR  RATIO  FOR  ISOHYET  A  AT  20000  SO.  MI.  “0.I06 


DEPTH  VS.  DURATION 


ISOHYET 

SMIN 

1.0MIM 

15MIN 

30MIN 

1-HR 

2 -HR 

3 -HR 

6-HR 

12-HR 

18-HR 

2  4 -HR 

30-HR 

36-HR 

42-HR 

48-HR 

54-HR 

60-HR 

66-HR 

72-HR 

A 

O.SI 

1.05 

1.57 

3.02 

5.32 

8.76 

11.94 

18.51 

22.97 

25.44 

27.15 

28.50 

29.61 

30,56 

31.39 

32.12 

32.77 

33.37 

33.91 

B 

0.49 

0.98 

1.46 

2.81 

4.96 

8.18 

11.16 

17.34 

21.68 

24.12 

25.83 

27.18 

28.29 

29,24 

30,07 

30,80 

31.45 

32.05 

32.59 

C 

0.45 

0.90 

1.34 

2.59 

4.58 

7.56 

10.32 

16.09 

20.29 

22.71 

24.42 

25.77 

26.88 

27.83 

28.65 

29.38 

30.04 

30.63 

31.17 

0 

0.42 

0.83 

1.24 

2.39 

4.22 

6.99 

9.54 

14.92 

19.03 

21.43 

23.14 

24.49 

25.60 

26,55 

27.37 

28.10 

28.76 

29.35 

29.89 

E 

0.18 

0.76 

1.13 

2.1.9 

3.87 

6.42 

8.76 

13.76 

17.79 

20.17 

21.88 

23.23 

24.34 

25.29 

26.11 

26.84 

27.50 

28.09 

28.64 

F 

0.15 

0.70 

1.05 

2.02 

3.58 

5.94 

8.11 

12.79 

16.75 

19.12 

20.83 

22.18 

23.29 

24.24 

25.06 

25,79 

26.45 

27.04 

27.58 

G 

o.n 

0.65 

0.96 

1.86 

3.29 

5.47 

7.46 

11.82 

15.70 

18.07 

19.78 

21.12 

22.24 

23.18 

24.01 

24.74 

25.39 

25.99 

26.53 

H 

0.30 

0.59 

0.88 

1.69 

3.00 

4.99 

6.82 

10.85 

14.66 

17.01 

18.72 

20.07 

21.18 

22.13 

22.96 

23,69 

24.34 

24.94 

25.48 

I 

0.27 

0.53 

0.79 

1.53 

2.71 

4.52 

6.17 

9.88 

13.64 

15.98 

17.69 

19.04 

20.15 

21.10 

21.92 

22,65 

23.31 

23,90 

24.45 

J 

0.24 

0.47 

0.70 

1.36 

2.42 

4.04 

5.52 

8.91 

12.61 

14.95 

16.66 

18.01 

19.12 

20.07 

20.89 

21,62 

22.28 

22.87 

23.41 

K 

0.21 

0.42 

0.62 

1.20 

2.15 

3.61 

4.94 

8.04 

11.66 

14.00 

15.71 

17.05 

18.17 

19.11 

19.94 

20.67 

21.32 

21.92 

22.46 

L 

0.19 

0.37 

0.55 

1.05 

1.89 

3.18 

4.36 

7.17 

10.75 

13.08 

14.79 

16.14 

17.25 

18.20 

19.02 

19.75 

20.41 

21.00 

21.54 

M 

0.07 

0.14 

0.2). 

0.41 

0.82 

1.61 

2.36 

4.26 

6.77 

8.54 

9.84 

10.87 

11.71 

12.43 

13.06 

13,61 

14.11 

14.56 

14,98 

N 

0.04 

0.08 

O.IV 

0.23 

0.45 

0.89 

1.31 

2.42 

4.10 

5.44 

6.41 

7.18 

7.81 

8.35 

8.82 

9.24 

9.61 

9.95 

10.26 

0 

0.02 

0.04 

0.06 

0.11 

0.22 

0.44 

0.65 

1.26 

2.32 

3.20 

3.84 

4.35 

4.77 

5.12 

5.43 

5.70 

5.95 

6.17 

6.38 

P 

o.ni 

O.Ol 

0.02 

0.03 

0.07 

0.13 

0.20 

0.39 

0.72 

1.00 

1.21 

1.37 

1.50 

1,62 

1.72 

1.81 

1.88 

1.96 

2.02 

Q 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

R 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

S 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

AVERAOE 

0.21 

0.42 

0.63 

1.21 

2.18 

3.69 

5.07 

8.24 

11.64 

13.84 

15.43 

16.69 

17.73 

18.62 

19.39 

20.07 

20.68 

21.24 

21.75 

4/ 


TAflLB  12(f) 

Sanole  Outout  from  HMR52  (contlnue'1) 


PROSABr.E  HWCrmiM  STORM  FOR  LHOM 


DAY  1 


TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

TIME 

PRBCTPTT&TTnN 

INCH 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

0200 

0.17 

0.17 

0800 

0.20 

0.71 

1400 

0.26 

1.38 

2000 

0.36 

2,24 

0400 

0.17 

0.34 

1000 

0.20 

0.92 

1600 

0.26 

1.63 

2200 

0.36 

2.S8 

OfiOO 

0.17 

0.61 

1200 

0.20 

1.12 

1800 

0.26 

1,89 

2400 

0.36 

2.93 

6-HR 

TOTAL 

0,61 

0.61 

0.77 

1.04 

DAY  2 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

INCH 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

0200 

0.49 

3.42 

0800 

0,95 

6.48 

1400 

2.67 

10,61 

2000 

0.63 

17.00 

0400 

0.63 

3.95 

1000 

1,12 

6,69 

1600 

3.69 

14.19 

2200 

0.72 

17.72 

0600 

0,68 

4.63 

1200 

1.34 

7.94 

1800 

1.98 

16.17 

2400 

0.64 

18.36 

6-HR 

TOTAL 

1.60 

3.41 

8.24 

2.19 

DAY  ■) 

TIME 

precipitation 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

TIME 

PRECIPIfATTON 

INCH 

total 

INCR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

0200 

0.42 

18.78 

0800 

0.30 

19.92 

1400 

0.23 

20.74 

2000 

0.19 

21.38 

0400 

0.42 

19.20 

1000 

0.30 

20.21 

1600 

0.23 

20.96 

2200 

0.19 

21.86 

0600 

0.42 

19.62 

1200 

0.30 

20.61 

1800 

0.23 

21.19 

2400 

0.19 

21.75 

6 -HR 

TOTAL 

1.26 

0.89 

0.68 

0.66 
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8 . 3  Prec  Ua U on  F lie  for  Ra Onfall-Runoff  Model  Usage 

Hyetographs  may  be  written  to  a  file  for  use  In  a  rainfall-runoff  model 
such  as  HEC-1. 

8 . ? . 1  Punch  File 

The  option  of  writing  hyetographs  to  a  punch  or  card-image  file  Is 
available.  Each  hyetograph  Is  preceded  by  a  line  giving  the  subbasin  name  and 
time  Interval  (see  Table  13.)  The  punch  data  must  be  merged  with  Input  data 
for  the  rainfall-runoff  model. 

8.2.2  DSS  File 

For  computer  systems  where  the  HEC  Data  Storage  System  (HECOSS)  Is 
available,  hyetographs  may  be  transferred  to  HEC-1  through  a  DSS  file.  See 
Section  9  for  an  example. 


Table  13 

"PUNCH"  Output  file  from  HNR52 


LEW  T'l’^ER'/AL  =  \10  Mt'l 

PI  0.l'i9  0  ln9  0.104  •).104 

PI  0.T47  0.147  0.43‘i  0.523  0.534 

PI  1.980  0.325  0.720  0.544  1.420 

PI  0.223  0.228  0.228  0.135  0.185 


0.204  0.257 
0.950  1.115 
0.420  0.420 
0.135 


0.257  0.257  0.147 
1.141  2.571  1.685 
0.295  0.295  0.295 
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Section  9 


EXAHPIF  APPLICATION 


^  1  lnlroduc  L1_q.n 

Fig.  25  shows  the  watershed  above  Jones  Reservoir.  HHR5?  Is  used  to 
develop  PHS  hyetographs  for  the  four  subbasins.  HEC-1  Is  used  to  calculate 
and  route  the  PMF  through  Jones  Reservoir.  The  hyetographs  are  transferred 
from  HMR5?  to  HEC-1  using  the  HEC  Data  Storage  System  (OSS).  The  HMR5?  Input 
data  for  this  example  are  shown  In  Table  14. 

9 . 2  HMR 52  IMU t  and  Output 

9.2.1  PMS  for  Jo  t a  1  Ba  s  I  n 

Data  for  total  basin  are  shown  on  lines  2-32  of  lable  14.  The 
boundary  of  the  total  watershed  was  digitized  and  entered  on  BX  and  BY  cards. 
Depth -area-duration  data  from  HMR  No.  51  have  been  put  on  HP  cards.  The 
preferred  storm  orientation  Is  on  an  HO  card. 

Storm-area  size  and  orientation  are  shown  as  zeroes  on  the  SA  card,  so 
HHR5?  will  select  the  storm-area  size  and  orientation  which  produce  maximum 
precipitation  on  the  watershed.  Tables  15(a)  and  15(b)  show  average  6-hr 
Incremental  depths  for  various  storm-area  sizes  and  orientations  used  In  the 
selection  process.  By  default,  HMR52  will  maximize  the  basin-average 
precipitation  for  the  three  largest  6-hr  Increments. 


Figure  25.  Jones  Reservoir  Watershed 
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TABLI?  14 

HMRSZ  Inout  Jones  Reservoir 
RBC  PROBABLE  MAXIMUM  STORM  INPUT  DATA 


[ME 

ID, . 

. 1 

. 7  .  . 

_ 3. 

_ 5, 

. .  .  .7  . . 

. R.  . 

....10 

1 

ID 

HMR52  INPUT 

UATA 

FOR  SAMPLE 

PMF  CALCULATION 

BN 

TOTAL 

1 

HS 

1 

4 

8X 

27.9 

27.4 

27.4 

26.9 

26.8 

26.5 

26.6 

26.2 

25.9 

25.4 

5 

8X 

26.1 

23.6 

21.4 

22.6 

22.5 

21.6 

21.3 

20.3 

19.9 

19.7 

6 

BX 

18.4 

18.2 

18.4 

18.4 

18.2 

19.3 

19.8 

20.7 

20.8 

21.8 

7 

ax 

22.4 

22.4 

23.0 

23.3 

23.7 

24.1 

24.5 

24.8 

25.3 

8 

BX 

26.4 

26.6 

27.5 

27.7 

27.8 

28,5 

29.6 

10.8 

11.2 

31.7 

8 

BX 

12.2 

14.6 

14.6 

34.9 

35.7 

36.4 

37.4 

38.2 

39,3 

39,6 

10 

BX 

39.4 

19.6 

18.6 

39.0 

40.0 

40.7 

41.4 

41.8 

42,2 

41.9 

11 

BX 

42.2 

41.4 

41.0 

43.0 

41.6 

41.0 

40.6 

40,5 

40.5 

40.1 

12 

BX 

19.6 

19.1 

18.9 

38.4 

38.0 

37.6 

36.7 

15,6 

15.4 

31.7 

11 

ax 

11.6 

10.9 

10.9 

30.7 

30.1 

29.9 

29.2 

28.2 

14 

BY 

89.1 

89.8 

90.1 

90.1 

90.8 

91.4 

92.3 

92.2 

92.3 

92.0 

10 

BY 

92.0 

90.6 

90.1 

89.9 

89,6 

89.2 

89.3 

87.9 

87.9 

87.3 

16 

BY 

86.6 

86.1 

86.9 

85.4 

86.0 

85.0 

85.4 

84.6 

83.9 

83.2 

17 

BY 

83.2 

82.9 

82.4 

81.4 

80.7 

80.4 

80.4 

80.0 

80.0 

79,8 

18 

BY 

79.6 

79.2 

78.8 

78.1 

77.2 

76.6 

76.1 

76,3 

76.2 

75.5 

1  <l 

BY 

75.4 

76.1 

74.4 

74.0 

73.8 

73.8 

73.9 

73.9 

74.2 

74.7 

20 

BY 

75.6 

76.4 

77,3 

78.0 

78.4 

79.0 

79,6 

80.3 

81,1 

82.0 

21 

BY 

81.8 

86.0 

86.6 

86.9 

87.1 

87.4 

87.7 

88.2 

90.6 

91.3 

22 

BY 

91.8 

92.1 

92.2 

92.2 

92.1 

91.5 

91.4 

91.8 

92.0 

92.0 

21 

BY 

90.9 

90.6 

90.1 

89.7 

89.7 

89.8 

89.2 

89.0 

24 

HO 

201 

25 

HP 

10 

26.5 

29.2 

32.0 

35.2 

36.8 

26 

HP 

200 

17.2 

20.5 

23.2 

26.2 

27.5 

27 

HP 

1000 

12.3 

15.6 

18.5 

21.1 

22.2 

HP 

6000 

7.5 

10.8 

13.2 

16.1 

17.2 

29 

HP 

10000 

5.7 

9.0 

11. 0 

13.5 

14.9 

30 

HP 

20000 

4.1 

7.1 

9.1 

11.8 

11.0 

11 

SA 

0 

0 

12 

ST 

120 

.10 

31 

BN 

1 

34 

BX 

18.2 

19.1 

19.8 

20.7 

20.8 

21.2 

21.8 

22.4 

22.8 

23.5 

35 

BX 

24.1 

24.6 

25.0 

26.2 

26.5 

27.0 

27.4 

27.4 

27.5 

27.9 

36 

BX 

27.9 

27.4 

27.4 

26.9 

26.8 

26.5 

26.6 

26.2 

26.9 

25.4 

37 

BX 

26.1 

21.6 

23.4 

22.6 

22.6 

21.6 

21.1 

20.3 

19.9 

19.8 

18 

BX 

19.7 

18.4 

18.2 

18.4 

18.4 

18.2 

39 

BY 

84.9 

86.0 

86.3 

84.6 

83.9 

83.6 

83.2 

83.2 

81.3 

82.7 

40 

BY 

R2.7 

81.4 

84.2 

84.5 

85.1 

84.8 

85.6 

87.9 

88.1 

88.5 

41 

BY 

89.1 

89.8 

90.1 

90.1 

90.8 

91.4 

92.1 

92.2 

92.1 

92.0 

42 

BY 

92.0 

90.6 

90.3 

89.9 

89.5 

89.2 

89.1 

87.9 

87.9 

87.6 

43 

BY 

87,1 

36.6 

86.1 

85.9 

95.4 

85.0 

44 

ZW 

A=SAMPLR 

4S 

m 

4r, 

BX 

27.4 

27.4 

77.0 

27.3 

27.4 

21,1 

27.9 

27.8 

27.2 

27.3 

47 

BX 

27.6 

2S.\ 

28.7 

29.6 

29.8 

11.0 

11.8 

12.1 

12.7 

13.0 

4R 

BX 

S 

13.9 

14.7 

15.2 

35.5 

16.6 

15.8 

16.1 

16.4 

37.0 
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TABLR  14 

HMR52  Inoufc  for  .Ton<?s  Reservoir  (continue'!) 


HBC  PROBABLR  MAXIMUM  STORM  INPtJT  DATA 


LI^JE 

ID. , 

. 7.. 

. 3. . 

. 4. . 

c; 

. 7.. 

. B. . 

. . . .10 

49 

BX 

. . 

38.1 

38.2 

38.5 

38.9 

39.4 

39.5 

19.7 

39.8 

40.2 

40.5 

sn 

BX 

40.5 

40.5 

40.5 

40.1 

39.5 

39.1 

38.9 

38.4 

38.0 

37.6 

51 

8X 

36.6 

35.6 

35.6 

35.4 

14.8 

13.7 

33.6 

33.5 

32.5 

30.9 

52 

BX 

30.9 

30.7 

30.1 

29.9 

29.2 

28.2 

27.9 

27.9 

27.5 

53 

BY 

87.9 

85.6 

84.8 

84.2 

83.6 

83.6 

82.6 

82.1 

81.1 

81,0 

54 

BY 

81.0 

80.9 

80.7 

80.7 

80.6 

80.7 

81.1 

81.7 

82.2 

82,5 

55 

BY 

82.4 

82.9 

83.3 

83.1 

83.6 

84.2 

84.8 

85.1 

85.2 

85,2 

56 

BY 

85.9 

86.4 

86.6 

86.6 

86.6 

86.5 

86.9 

87.8 

88.2 

88.2 

57 

BY 

89.0 

89.8 

90.6 

91.3 

91.8 

92.1 

92.2 

92.2 

92,1 

91.5 

58 

BY 

91.5 

91,4 

91,8 

92.0 

92.0 

91.9 

91.1 

90.9 

90.8 

90.6 

59 

BY 

90.1 

89.7 

89.7 

89.8 

89,2 

89,1 

89.0 

89.1 

88.5 

88.1 

60 

ZW 

61 

BH 

3 

62 

BX 

27.6 

27.1 

27.6 

27,7 

27.7 

27.5 

27.7 

27.8 

28.5 

29.6 

63 

BX 

30.3 

30,8 

31,2 

11,7 

32.2 

14.5 

34.6 

14.9 

35.7 

16.4 

64 

BX 

37.4 

37.9 

18.2 

39.1 

39.1 

19.6 

39.4 

19.5 

39.1 

18.6 

65 

BX 

38.8 

39. n 

40.0 

40.7 

41.4 

41.8 

42.2 

41.9 

42.0 

42.2 

66 

BX 

42.8 

43.1 

43.0 

41.0 

42.4 

41.6 

41.0 

40.6 

40.5 

40.1 

67 

BX 

39.7 

39.7 

19.5 

39,0 

38.4 

18.2 

38.1 

17.0 

36.2 

35.9 

68 

BX 

35.7 

35.5 

35.1 

34.8 

31.9 

33.5 

33.0 

32.5 

12.3 

32.1 

69 

BX 

31.7 

11.2 

10,8 

29.5 

28,4 

28.1 

27.6 

70 

BY 

80.9 

80,1 

79.8 

79.5 

78.8 

78.8 

78.1 

77.2 

76.6 

76.1 

71 

BY 

76.2 

76.1 

76.2 

75.5 

75.4 

75.3 

74.4 

74.0 

73.8 

73.8 

72 

BY 

73.9 

73.9 

73.9 

74,1 

74.2 

74.7 

75.6 

76.3 

76.8 

77.3 

73 

BY 

77.6 

78.0 

78.4 

79,0 

79.6 

80,1 

81.1 

81.9 

83.1 

83.8 

74 

BY 

84,3 

85.0 

86.5 

86.9 

87.0 

87.1 

87.4 

87.7 

88.2 

88.1 

75 

BY 

87.7 

87.0 

86.5 

86.7 

86.5 

86.4 

85.9 

85.2 

85.1 

84.9 

76 

BY 

84.5 

83.6 

83.3 

83.3 

82.9 

82.5 

82.5 

82.1 

81.9 

81.5 

77 

BY 

81.1 

80.8 

80.7 

80.6 

80.7 

80.8 

80.9 

78 

zw 

79 

BN 

4 

80 

BX 

27.6 

27.1 

27.6 

27.7 

27.7 

27.5 

26.9 

26.5 

26.0 

25.4 

81 

BX 

25.3 

24.8 

24.5 

24.1 

21.7 

23.1 

23.0 

22.7 

22.1 

22.2 

82 

BX 

22.4 

22.8 

23.5 

24.0 

24.6 

25.0 

25.2 

26.5 

27.0 

27.3 

83 

BX 

27.4 

27.7 

27.9 

27.8 

27.2 

27.3 

27.6 

84 

BY 

80.9 

80.1 

79.8 

79.5 

78.8 

78.8 

79.0 

79.2 

79.3 

79.5 

85 

BY 

79.8 

80.0 

80.0 

80.4 

80.4 

80.7 

81.4 

81.9 

82.3 

82.9 

86 

BY 

83.2 

81.1 

82.7 

82.7 

83.4 

84.2 

84.5 

85.1 

84.2 

87 

BY 

83.6 

83.6 

82.6 

82.1 

81.3 

81.0 

81.0 

88 

ZW 

89 

zz 
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In  Table  15(a),  storm  precipitation  Is  calculated  for  various  storm-area 
sl/es.  Storm  orientation  was  chosen  for  this  table  by  minimizing  the  moment  of 
Inertia  about  the  major  axis  of  the  elliptical  storm  pattern.  In  Table  15(b) 
the  storm-area  size  Is  fixed  at  the  size  which  produced  maximum  precipitation 
on  the  watershed  In  Table  15(a),  and  orientation  Is  varied  from  140  to  310  by 
ten-degree  Increments.  The  last  two  storms  In  lable  15(b)  have  orientations 
which  are  five  degrees  to  either  side  of  the  best  previous  orientation.  By 
coincidence,  the  best  orientation  Is  285  ,  so  the  last  two  lines  are  repeats 
of  storms  calculated  previously. 

9.2.2  Subba sj^n  Rye tographs 

Data  for  the  four  subbasins  follows  the  data  for  the  total  watershed 
In  the  HMR5?  Input,  Table  14.  HMR52  uses  the  same  storm  specifications  as 
were  developed  for  the  total  basin  to  calculate  portions  of  the  storm 
occurring  over  each  subbasin.  These  storms  are  shown  In  Tables  16  through  19. 

A  ZW  card  Is  Included  with  data  for  each  subbasin  to  Indicate  that  the 
storm  hyetographs  are  to  be  written  to  a  DSS  file. 

NOIt  -  HtCOSS  has  been  developed  for  Corps  of  Engineers  use  at  selected 
computer  sites.  DSS  subroutines  are  not  generally  available  to  non-Corps 
users.  Users  without  access  to  OSS  should  use  the  punch  (tape)  file  option  of 
the  program,  see  section  8.2.1. 

9 . 3  NEC- L  Input  and  .0utp.ut 

The  rainfal  1 -runoff  model  HEC-1  was  used  to  calculate  the  runoff  from  the 
storm  hyetographs.  Table  20  shows  HFC-1  Input  data  for  this  application,  and 
Table  21  shows  HEC-1  summary  output  for  the  storm  generated  by  HMR5?. 

9.4  PHF  Q#icuJ.atJo_n 

Ihe  objective  of  calculating  a  PMf  Is  to  obtain  the  largest  flood  which 
may  reasonably  occur  over  the  watershed.  Because  of  hydrologic 
characteristics  of  a  watershed,  the  largest  flood  may  not  result  from  the 
storm  which  produces  the  most  precipitation. 

Several  trials  were  made  to  calculate  the  PMf  for  Jones  Reservoir. 

Results  from  these  trials  are  tabulated  In  lable  22. 

Trial  I:  used  storm-area  size  and  orientation  selected  by  HMR52;  default 
value  for  position  of  peak  6  hr  Interval;  basin  centroid  for  storm  center 
location,  see  fig.  26. 

Trial  2:  same  as  trial  1,  except  peak  6  hr  Interval  Is  shifted  to  the 
10th  position  (hours  54-60).  This  raised  the  peak  flow  slightly  so  the  10th 
position  was  used  for  subsequent  trials,  except  trial  6. 

Trial  3:  same  as  trial  2,  except  the  Isohyelal  pattern  was  centered  on 
the  watershed  by  visual  adjustment  of  the  storm  center. 

Irlal  4;  same  as  trial  3,  except  a  smaller  storm  area  sl/e  was  used. 
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Trial  5;  Since  subbasins  1,  ?,  and  3  produce  most  of  the  runoff,  the 
storm  was  centered  over  these  subbasins. 

Trial  same  as  trial  1,  except  the  order  of  the  6-hr  Increments  Is  as 
shown  In  Fig.  7. 

9 . 5  Conclusions  and  Recommendations 

The  adopted  PMF  for  this  watershed  would  be  the  flood  calculated  In  trial 
?.  The  results  from  trial  1,  using  program  defaults,  could  easily  be  adopted 
for  the  PMF,  since  the  difference  Is  less  than  0.4  percent  for  peak  Inflow  and 
less  than  0.7  percent  for  peak  outflow. 

In  most  cases  default  values  In  HMRS2  will  suffice  to  determine  the 
PMF.  Exceptions  will  occur  where  the  basin  has  an  Irregular  shape  or  unusual 
hydrologic  characteristics.  Consideration  of  the  following  characteristics 
may  be  useful  In  selecting  parameters  for  alternative  trials: 

relative  timing  of  runoff  from  subbasins 

direction  of  major  stream 

proximity  of  storm  center  to  a  major  stream 

relationship  between  storm-area  Isohyet  and  basin  boundary 


Figure  ?6.  Storm  Pattern  For  trials  1.  ?.  and  b 
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TABLE  16 

Probable  Maximum  Storm  for  Subbasin  1 


DAY  1 


time 

PRECIPITATION 

TIME 

PRECIPITATION 

TIME 

INCR 

TOTAL 

INCR 

total 

0200 

O.ll 

0.11 

0800 

0.13 

0.46 

1400 

0400 

O.ll 

0.22 

1000 

0.13 

0.59 

1600 

0600 

o.n 

0.33 

1200 

0.13 

0.72 

1800 

6 -HR 

TOTAL 

0.13 

0.40 

DAY  2 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

TIME 

INCR 

TOTAL 

INCR 

total 

0200 

0.34 

2.28 

0800 

0.82 

3.91 

1400 

0400 

0.38 

2.66 

1000 

1.05 

4.95 

1600 

0600 

0.41 

1.09 

1200 

1.37 

6.32 

1800 

6 -HR 

TOTAL 

1.1.5 

1.24 

DAY  1 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

TIME 

INCR 

total 

INCR 

total 

0200 

0.29 

23.82 

0800 

0.20 

24.60 

1400 

0400 

0.29 

24.11 

1000 

0.20 

24.79 

1600 

0600 

0.29 

24.40 

1200 

0.20 

24.99 

1800 

PRECIPITATION 

TIME 

PRECIPITATION 

INCR 

TOTAL 

INCR 

TOTAL 

0.17 

0.89 

2000 

0.23 

1.46 

0.17 

1.06 

2200 

0.23 

1.70 

0.17 

1.23 

2400 

0.23 

1.93 

0.51 

0.70 

PRECIPITATION 

TIME 

PRECIPITATION 

INCR 

TOTAL 

INCR 

TOTAL 

3.94 

10.27 

2000 

0.67 

22.50 

8.99 

19.26 

2200 

0.55 

23.05 

2.57 

21.83 

2400 

0.48 

23.53 

15.51 

1.69 

PRECIPITATION 

TIME 

PRECIPITATION 

INCR 

TOTAL 

INCR 

total 

0.15 

25.14 

2000 

0.12 

25.55 

0.15 

25.29 

2200 

0.12 

25.67 

0.15 

25.44 

2400 

0.12 

25.79 

6-HR  total  D.S7 


0.59 


0.45 


0.36 


TABLE  17 

Probable  Maximum  Storm  for  Subbasin  2 


TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

INCH 

TOTAL 

INCR 

total 

INCR 

TOTAL 

INCR 

TOTAL 

0200 

0.11 

0.11 

0800 

0.13 

0.44 

1400 

0.16 

0.87 

2000 

0.21 

1.42 

0400 

O.ll 

0.21 

1000 

0.13 

0.67 

1600 

0.16 

1.03 

2200 

0.23 

1.65 

0600 

O.ll 

0.12 

1200 

0.13 

0.70 

1800 

0.16 

1.19 

2400 

0.21 

1.88 

6-HP 

TOTAL 

0.12 

0.19 

0.49 

0.68 

DAY  2 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

INCR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

0200 

0.11 

2.21 

0800 

0.81 

1.81 

1400 

4.13 

10.36 

2000 

0.66 

22.64 

0400 

0.37 

2.68 

1000 

1.04 

4.86 

1600 

8.94 

19.29 

2200 

0.64 

23.18 

0600 

0.42 

1,00 

1200 

1.18 

6.21 

1800 

2.69 

21.98 

2400 

0.47 

23.64 

6-HR 

TOTAL 

1.12 

3.23 

16.76 

1.66 

DAY  ^ 

TIME 

PRECIPITATION 

TIME 

precipitation 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

INCR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

0200 

0.2S 

23.93 

0800 

0.19 

24.68 

1400 

0.14 

26.21 

2000 

0.12 

25.61 

0400 

0.28 

24.21 

1000 

0.19 

24,87 

1600 

0.14 

25.16 

2200 

0.12 

25.73 

0600 

0.28 

24.49 

1200 

0.19 

26.06 

1800 

0,14 

26.50 

2400 

0  12 

26.84 

6-HR 

TOTAL 

0.8S 

0.67 

0.43 

0.15 

'jI 


TABLE  18 

Probable  Maximum  Storm  for  Suhhasln  3 


DAY  1 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

TIME 

precipitation 

TIME 

PRECIPITATION 

INCR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

0200 

O.ll 

O.ll 

0800 

0.13 

0.4S 

1400 

0.16 

0.87 

2000 

0.23 

1.43 

0400 

o.ll 

0.21 

1  000 

0.13 

O.SR 

1600 

0.15 

1.03 

2200 

0.23 

1.66 

0800 

O.ll 

0.32 

1200 

0.13 

0.70 

1800 

0.1.6 

1.20 

2400 

0.23 

1,89 

fi-HR 

TOTAL 

0.32 

0.39 

0.49 

0.69 

DAY  2 

TIME 

PRECIPITATION 

TIME 

PRECIPC'T’ATTON 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

TNCR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

0200 

0.34 

2.22 

0800 

0.80 

3.81 

1400 

3.93 

10.10 

2000 

0.65 

21.77 

0400 

0.37 

2.69 

1000 

1.02 

4.83 

1600 

8.44 

18.64 

2200 

0.54 

22.30 

0600 

0.42 

3,01 

1  200 

1.34 

6.17 

1800 

2.67 

21.11 

2400 

0.47 

22.77 

«'HR 

TOTAL 

1.12 

3.16 

14.94 

1,65 

DAY  3 

TIME 

PRECIPITATION 

TIME 

precipitation 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

INCR 

total 

INCR 

TOTAL 

INCR 

total 

INCR 

TOTAL 

0200 

0.28 

23.06 

0800 

0.19 

23.81 

1400 

0,14 

24.34 

2000 

0,12 

24.75 

0400 

0.28 

23.34 

1  000 

0.19 

24.00 

1600 

0.14 

24.48 

2200 

0.12 

24,86 

0600 

0.28 

23.62 

1200 

0.19 

24.70 

1800 

0.14 

24.63 

2400 

0,12 

24.98 

8-HR 

TOTAI, 

0.86 

0.67 

0.43 

0.35 

!)0 


TABLE  19 

Prohahl<»  Maximum  Storm  for  Subbasin  4 


DAY  I 


TTMR 

PRBCtPITATICW 

TIME 

PRBCIPITATIOM 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

rvrR 

TOTAL 

I  NCR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

0200 

O.ll 

0.11 

0800 

0.11 

0.46 

1400 

0.17 

0.89 

2000 

0.24 

1.47 

0400 

0.11 

0.22 

1000 

0.11 

0.69 

1600 

0.17 

1.06 

2200 

0.24 

1.70 

o#^no 

0.11 

0.11 

1200 

0.11 

0.72 

1800 

0.17 

1,23 

2400 

0.24 

1.94 

6 -HR 

TOTAL 

0.13 

0.40 

0.61 

0.71 

DAY  7 

TTME 

PRfriPTTATION 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

IHCR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

0200 

O.IS 

2.28 

0800 

0.82 

1.92 

1400 

3.96 

10.10 

2000 

0.67 

22.63 

0400 

O.IR 

2.66 

1000 

1.06 

4.97 

1600 

9.08 

19.18 

2200 

0.66 

23.18 

O^iOO 

0.41 

1.09 

'  200 

1.18 

6.14 

1800 

2.68 

21.96 

2400 

0.48 

21.66 

6 -HR 

TOTAL 

1.16 

1.26 

16.62 

1.70 

DAY  1 

TIME 

PRVy'TPtTA'^TON 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

TIME 

PRECIPITATION 

iNrR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

INCR 

TOTAL 

0200 

0.29 

21.96 

0800 

0.20 

24.71 

1400 

0.16 

26.27 

2000 

0.12 

26.69 

0400 

0.29 

24.24 

1000 

0.20 

24.91 

1600 

0.16 

26.42 

2200 

0.12 

25,81 

OfiOO 

0.29 

24.61 

1200 

0.20 

26.12 

1800 

0.16 

2.6.67 

2400 

0.12 

26.91 

fi-HR  TOTAl,  n.RR 


0.S9 


0.4S 


0.16 


LINE 

1 

2 

3 

4 

6 

7 

8 
9 

10 

11 

12 

13 

14 
18 
16 

17 

18 

19 

20 

21 

22 

23 

24 
28 
26 

27 

28 
29 
10 

31 

32 

33 

34 
38 

36 

37 

38 

39 
49 
41 
4? 

43 

44 
48 
48 


HEC-l  Input  for  -Tones  Reservoir  Watershe'3 


HEC-1  INPUT 

ID . 1 . 2 . 3 4 8 6 7 8 9 . 10 

ID  HRCl  DMA  FOR  SAMPLE  PMF  CALCULATION 

IT  120  01.TAN99  0000  80 

10  8 

EK  1 

KM  BEAR  CREEK  IMFI/DM  INTO  JONES  RBSERVOUR 

BA  81 . 80 

ZR  =PI  A’SAMPLE  C=PRECIP-INC 
LU  1  .06 

UR  9.00  .48 

BF  -3.  -.35  1.0092 

RK  4 

KM  LOCAL  RUNOFF  INTO  JONES  RESERVOIR 

BA  20.30 

ZR  =PI 


WI 

US 

UA 

1 

2.60 

10 

.06 

.48 

13 

10 

20 

29 

40 

81 

70 

82 

93 

100 

KK 

KM 

1+4 

BEAR 

CREEK  AND 

LOCAL 

TMFLOM 

COMBINED 

HC  2 

KK  2 

KM  DEER  CREEK  INTO  JONES  RESERVOIR 
BA  98.60 

ZR  =P1 

LU  1  . 06 

OS  11.00  .45 

>'A 

KK  3 

KM  SKUNK  HOLLOW  INTO  JONES  RES 
BA  116.90 

ZR  =P1 

US  14.00  .48 

KK  2-t3 

KM  COMBINED  FLOW  FROM  DEER  CR  AND  SKUNK  HOLLOW 

HC  2 

KK  INFLOW 

KM  TOTAL  INFLOW  INTO  JONES  RES 
HC  2 

KK  JONES 

KM  ROUTE  THROUOH  JONES  RES 


PS 

1 

PLOW 

-1 

S7 

1  4*^000 

iinoo 

14^^00 

22000 

36000 

54000 

82000 

110700 

118000 

125000 

n4ooo 

S'? 

sn 

0 

i?sono 

17000 

isoooo 

1  3000 

14000 

isnoo 

16000 

17000 

30000 

54000 

flfiono 

7.1 


pr>.m  1 
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Table  21 


HEC-1  Summary  Output  for  Trial  1 


RUNOFP  SUMMARY 

PUM  IN  CUBIC  FEBT  PER  SECOND 
TIME  IN  HOURS,  AREA  IN  SQUARE  MILES 


OPERATION 

STATION 

PEAK 

FLOW 

TIME  OF 
PEAK 

AVERAGE 
6 -HOUR 

FLOW  FOR  MAXIMUM 
24-HOtJR 

PERIOD 

72-HOOR 

BASIN 

AREA 

HYDROGRAPH 

AT 

1 

29528. 

48.00 

28056. 

20818. 

11805. 

51.80 

HYDROGRAPH 

AT 

4 

26798. 

42.00 

21758. 

12922. 

8339. 

20.30 

2  COMBINED 

AT 

1+4 

43374. 

44.00 

42278. 

33231. 

20144. 

72.10 

HYDROGRAPH 

AT 

2 

48055. 

50.00 

46176. 

36320. 

20357. 

98.60 

HYDROGRAPH 

AT 

3 

44777. 

52.00 

43877. 

36318. 

20591. 

116.90 

2  COMBINED 

AT 

2+3 

90924. 

52.00 

88749. 

72386. 

40948. 

215.50 

2  COMBINED 

AT 

INFLOW 

127493. 

50.00 

124553. 

103675. 

61092. 

287.60 

ROUTED  TO 

.TONES 

94664. 

60.00 

91656. 

74403. 

42765. 

287.60 

Table  22 

Summary  of  PMf  Calculations 


Trial 

Position  of 
Peak  6-hr 

Interval 

Storm  Area 
ml  2 

Orientation 

degrees 

Storm 

X 

miles 

Center 

y 

miles 

Total 
Ralnfal 1 

Inches 

Peak 

Inflow 

cfs 

1 

7 

300 

285 

32.2 

83.8 

25.50 

127,500 

2 

10 

300 

285 

32.2 

83.8 

25.50 

128,000 

3 

10 

300 

285 

31.0 

83.6 

25.44 

127,650 

4 

10 

175 

290 

31.0 

83.6 

24.86 

125,200 

5 

10 

300 

296 

32.7 

84.0 

25.41 

127,200 

6 

9* 

300 

285 

32.2 

83.8 

25.50 

125,100 

♦Temporal  distribution  based  on  Fig.  7. 


Peak 

Outflow 

cf  s 

94.650 

95.300 
95,250 

91.650 
93,900 

94.300 
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Section  10 


COMPUTER  REQUIREMENTS 


10. 1  Source  Language 

HMR52  source  code  1s  written  In  FORTRAN  77.  It  must  be  compiled  with  a 
Fortran  77  compiler  or  a  Fortran  compiler  which  can  Interpret  IF-THEN-ELSE- 
ENDIF  statements. 

10.2  Core  Storage 

Program  HMR52  uses  220K  bytes  on  a  Harris  500  computer  or  45K  words 
(decimal)  on  a  COC  7600.  Table  23  gives  storage  requirements  and  compile  and 
execution  times.  Data  for  other  computers  will  be  added  to  the  table  as  It  Is 
supplied  by  users. 


Table  23 


Memory  and  Time  Requirements 


Computer 

Harris  500 

CDC  Cyber 

Memory 

220K  bytes 

45K  words 

Compile  time 

93 

6.57 

(CPU  sec) 

Execution  time* 

43 

16.3 

(CPU  sec) 

•Execution  time  Is  for  example  data  In  Section  9,  Table  14. 
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HMR5?  uses  6  disk  files: 


Input,  card  or  card  Image; 

output,  formatted  for  132  column  printer; 

punch; 

two  scratch  files;  and 
a  table  file. 

File  characteristics  are  listed  In  Table  24. 

The  table  file  (unit  INP52)  contains  data  listed  In  Tables  1-5.  Data  on 
this  file  are  read  Into  the  program's  working  storage  at  the  beginning  of  each 
execution. 

10.4  HEC-Supplled  Magnetic  Tape 

NEC  programs  are  distributed  on  half  Inch  by  2400  foot  magnetic  tapes. 
Files  on  the  tape  for  HMR52  are  the  source  code,  a  file  containing  data  from 
Tables  1-4  and  10,  sample  Input  data,  and  the  current  Input  description  for 
HMR52. 

Ihe  source  code  on  the  tape  conforms  to  ANSI  standards  for  Fortran  77. 

Ihe  code  Is  not  Intended  for  a  particular  computer,  although  minor  changes  may 
be  necessary  to  use  the  program  on  different  machines. 

Ihe  sample  Input  data  are  the  same  as  used  In  Section  9.  This  data  may 
be  used  to  verify  that  the  program  Is  working  correctly. 

10.5  Machine-Specific  Code 

The  only  code  In  the  HMR52  source  code  which  Is  specific  to  a  particular 
computer  Is  the  use  of  subroutines  DAlt  and  TIME  In  subroutine  BANNER.  These 
subroutines  return  the  current  date  and  time.  This  Information  Is  printed  on 
the  banner  output  page  to  Identify  when  a  run  Is  made.  Some  code 
modifications  may  be  required  since  subroutines  DAlt  and  TIME  are  not  standard 
on  all  machines.  These  subroutine  calls  can  be  removed  without  affecting  the 
program's  performance. 


Table  24 


File  Characteristics 


Fortran 
Unit  No. 

Variable 

Name 

Description 

Maximum 

Record  Length 

5 

INP 

Input  data 

80  characters 

6 

IPRT 

Printer  output 

132  characters 

7 

IPNCH 

File  output* 

80  characters 

19 

SCRTCH 

scratch;  used  to  convert 
time  from  411  to  A4  format 

5  characters 

20 

1C 

working  data  file;  contains 
reformatted  Input  data  with 
line  number  and  next  card 

Is  appended  to  beginning 
of  each  card 

89  characters 

52 

INP52 

contains  data  tables  from 

HHR  No.  52 

80  characters 

*  This  file  maybe  punched  or  saved  on  disk/tape  for  use  as  Input  to  a 
rainfall-runoff  model. 
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APPENDIX  A 


HMR52  INPUT  DESCRIPTION 
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BA  CARD  -  Basin  Area . 70 
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GO  CARO  -  Compute  Storm  with  Current  Data . 72 

HO  CARD  -  Preferred  Orientation . 73 

HP  CARD  -  Depth-Area-Ouratlon  Data . 73 

ID  CARD  -  Job  Identification . 74 

PL  CARD  -  Plot  Control . 75 

PU  CARD  -  Punch  PMS . 76 

SA  CARD  -  Storm  Area  and  Orientation . 77 

SC  CARD  -  Storm  Center  Coordinates . 78 

SD  CARD  -  Arrangement  of  6-hr  Increments  In  In  PHS . 78 

Temporal  Distribution 

SI  CARD  -  Storm  Isohyet  Areas  (two  cards) . 79 

ST  CARD  -  Temporal  Distribution . 80 

ZH  CARD  -  Write  PMS  to  DSS . 81 

ZZ  CARO  -  End  of  Data . 8? 

Example . 83 
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HMR6?  INPUT  DESCRIPIIOK 


INPUT  FORMAT 

Input  data  Is  read  from  cards  or  card  Images.  Each  card  Is  divided  Into  a 
two-coTumn  card  Identification  field  (columns  T  and  ?)  and  ten  data  fields.  The 
first  data  field  has  six  columns  (columns  3-8).  The  remaining  nine  fields  each 
have  eight  columns. 

field  0  (zero)  Is  the  card  Identification  field  (columns  1  and  ?).  Field  1 
Is  the  first  data  field  (beginning  In  column  3). 

Under  the  value  heading  a  plus  (♦)  Is  used  to  Indicate  where  numeric  values 
should  be  entered.  A  blank  numeric  value  Is  Interpreted  as  a  zero.  (AN) 
Indicates  that  alphanumeric  characters  may  be  entered. 

Data  may  be  entered  on  the  cards  In  either  FREE  FORMAl  or  FIXED  FORMAT 
(default).  The  FORMAT  (fixed  or  free)  Is  controlled  by  a  switch  which  can  be 
changed  at  anytime  by  inserting  a  *FIX  or  a  ‘FREE  card  In  the  data  deck. 

For  FREE  FORMAT  each  Item  Is  separated  by  one  or  more  blanks  or  a  comma.  A 
blank  field  Is  designated  by  two  successive  commas.  There  may  be  more  or  less 
than  ten  coordinates  on  a  BX  or  BY  card.  The  FREE  FORMAl  reader  will  treat 
these  cards  as  If  they  were  on  continuous  cards. 

for  FIXED  FORMAT  each  number  which  does  not  contain  a  decimal  point  must  be 
right  justified  In  Its  field. 

A  ‘FIX  card  Indicates  that  the  following  data  cards  use  FIXED  FORMAT,  and 
a  ‘FREE  Indicates  FREE  FORMAT.  The  format  may  be  changed  as  often  as  desired 
and  at  any  location  In  the  data  deck. 

Comments  may  be  placed  anywhere  In  the  Input  deck  on  cards  with  double 
asterisks  (“)  In  the  card  Identification  field. 
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HHR52  INPUT  DESCRIPTION 


CARD  SEQUENCE 

Cards  are  grouped  by  drainage  basin.  A  BN  card  Is  used  to  Indicate  the 
beginning  of  data  for  a  basin.  The  next  cards  are  the  BS,  BX,  and  BY  cards 
which  describe  the  basin  boundary.  Next  come  the  HO  and  HP  cards  which  have 
meteorological  data  from  HMR  No.  52  and  HHR  No.  51,  then  the  SA,  SC,  SD,  and  ST 
cards  containing  data  for  a  particular  storm. 

The  cards  for  a  drainage  basin  may  be  followed  by  another  group,  beginning 
with  a  BN  card  for  a  different  drainage  basin,  or  a  GO  card  followed  by 
different  storm  data  for  the  current  basin,  or  a  ZZ  card  Indicating  the  end  of 
data  for  this  Job. 

Within  a  card  group  (beginning  with  BN  or  GO  card)  the  cards  may  be  placed 
In  any  order.  However,  It  Is  helpful  to  use  a  consistent  order  such  as 
alphabetical  order. 


DATA  REPETITION 


Once  data  has  been  read  by  the  program  It  Is  retained  In  memory  until  a  new 
card  of  the  same  type  Is  read.  For  example,  the  HO  or  HP  cards  need  only  occur 
In  the  data  for  the  first  drainage  basin  and  not  repeated  for  each  subsequent 
basin.  Similarly,  once  the  storm  area  and  orientation  have  been  established, 
they  will  be  used  for  all  subsequent  calculations  until  a  new  SA  card  Is  read. 

The  only  cards  which  must  be  Included  for  each  new  drainage  basin  are  the 
BN,  BX,  and  BY  cards.  All  other  cards  may  be  omitted  for  a  new  basin  and  the 
storm  will  be  calculated  using  data  from  the  previous  basin. 
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HMR52  INPUT  DESCRIPTION 


BA 

BA  CARD  -  BASIN  AREA 

This  card  Is  an  optional  card  which  Is  placed  after  the  BN  card.  A  BA  card 
Is  required  when  a  storm  Is  calculated  using  Isohyet  areas  on  SI  cards  and 
boundary  coordinates  are  not  given.  The  area  on  a  BA  card  Is  not  used  If 
boundary  coordinates  are  given. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTIDN 

0 

ID 

KA 

Card  Identification. 

1 

BARE  A 

t- 

Subbasin  area  In  square 

miles . 

BN  CARD 

-  BASIN  NAME 

The  BN  card  Is 
Identifies  the  basin 

placed 

at  the  beginning  of  data 

for  each 

drainage  basin 

FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

II) 

BN 

Card  Identification. 

1 

NAME 

AN 

Alphanumeric  name  for 

drainage 

basin  described 

data  on  following  cards.  (Maximum  8  characters  for 
FREE  FORMAT,  6  characters  for  FIXED  FORMAT.) 
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HMR52  INPUT  DESCRIPTION 

BS 

BS  CARD 

-  SCALE 

BX 

FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

BS 

Card  Identification. 

1 

SCALXY 

f 

Scale  of  boundary  coordinates 
coordinate  unit. 

In  miles  per 

BX  CARO  -  BOUNDARY  COORDINATES 

BX  and  BY  cards  contain  X  and  Y  coordinates  For  points  on  the  basin 
boundary.  The  points  may  be  entered  In  either  clockwise  or  counter-clockwise 
direction.  The  beginning  point  should  not  be  repeated  since  the  program 
automatically  closes  the  boundary. 

The  program  counts  the  number  of  boundary  points  by  finding  the  last 
non-zero  value  on  a  series  of  BX  or  BY  cards. 

IF  both  BA  and  SI  cards  are  used,  BX  and  BY  card  may  be  omitted. 

riELD  VARIABLE  VALUE  DESCRIPTION 

0  ID  BX  Card  Identification. 

X-coordInates  of  drainage  basin  boundary. 
Corresponding  to  y-coord1nates  on  BY  cards. 
(Maximum  -  TOO  values. ) 


1-10  XB 


HMRb2  INPUT  DESCRIPTION 


BY 

GO 


BY  CARO 

-  BOUNDARY 

COORDINATES 

FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

BY 

Card  Identification. 

1  -10 

YB 

> 

Y-coordInates  of  drainage  basin 

corresponding  to  x-coordinates  on 

(Maximum  -  TOO  values.) 


boundary 
BX  card. 


GO  CARD  -  COMPUTE  STORM  WITH  CURRENT  DATA 

The  GO  card  Is  used  to  Indicate  the  end  of  data  for  a  PMS  calculation.  The 
computer  will  calculate  the  PMS  using  the  current  data  before  reading  the  next 
card.  A  GO  card  Is  not  required  If  the  next  card  Is  a  BN  card  or  a  ZZ  card. 

FIELD  VARIABLE  VALUE  DESCRIPTION 

0  ID  GO  Card  Identification. 
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HMR52  INPUT  DESCRIPTION 


HO 

HP 


HO  CARD 

-  PREFERRED 

ORIENTATION 

FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

HO 

Card  Identification. 

1 

PORNT 

¥ 

Preferred  probable  maximum  storm  orientation  from 
HMR  No.  52. 

HP  CARD  -  DEPTH-AREA-OURATION  DATA 

Use  one  HP  card  for  each  storm  area  (10,  200,  1000,  5000,  10000,  and  20000 
square  miles).  Six  HP  cards  are  required,  one  for  each  area. 

FIELD  VARIABLE  VALUE  DESCRIPTION 

0  10  HP  Card  Identification. 

1  AREA51  *■  Storm  area  from  HMR  No.  51  for  PMP  In  Fields  2-6. 

2  H510AD(1,J)  *■  Probable  maximum  precipitation  for  6  hour  duration 

for  storm  area  1n  Field  1,  from  HMR  No.  51. 

3-6  H51DAD(1,J)  f  Similar  to  Field  2  for  durations  of  12,  24,  48  and 

72  hours. 
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HMR5?  INPUT  DESCRIPTION 


ID  CARD  -  JOB  IDENTIFICATION 

The  contents  of  this  card  are  read  and  printed  immediately  after  being 
read.  There  Is  no  limit  on  the  number  of  ID  cards.  ID  cards  are  usually  placed 
at  the  beginning  of  a  data  set  to  provide  title  Information.  They  may  be  placed 
within  a  data  set  to  describe  Individual  basins  or  storm  calculations. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

ID 

Card  Identification. 

1-10 

TITLE 

AN 

Alphanumeric  job  title  or  description 

PL  CARO  -  PLOT  CONIROL 


PL 


The  Pt  card  Is  used  to  control  printer  plots  of  the  drainage  basin  boundary 


and  PMS 

Isohyets . 

f  lEM 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

PL 

Card  Identification. 

1 

PLTYP 

0 

No  plots. 

1 

Plot  drainage  basin  boundary. 

2 

Plot  drainage  basin  boundary  with  PMS  Isohyets. 

3 

Make  two  plots,  a  plot  of  drainage  basin  boundary 
and  a  plot  of  the  boundary  with  PMS  Isohyets. 

2 

CHRPIN 

f 

Characters  per  Inch  for  printer  plot.  Default  Is 
10.) 

3 

LINPIN 

f 

Lines  per  Inch  for  printer  plot.  (Default  Is  6.) 

/S 


PU 


HNR52  INPUT  DESCRIPTION 


PU  CARO  -  SAVE  PNS  ON  FILE 


The  PU  Card  U  used  to  control  writing  the  PMS  to  a  file.  The  default 
format  Is  (2HPI,F6.2.9F8.2) . 


FIELD 

VARIABLE 

VALUE 

DESCRIPTICN 

0 

ID 

PU 

Card  Identification 

1 

PUNCH 

ON 

Write  PMS  precipitation  to  save  file 

OFF 

Do  not  write  PMS  to  save  file 
(default) . 

2-10 

PUNFMT 

AN 

Format  to  be  used  In  writing  to  save 
file.  If  Fields  2-10  are  blank  the 
previously  defined  format  will  be 
used. 

NOTE:  For  program  files  maintained  by  HEC  on  HARRIS  the  save  file  Is 
U8  by  default,  and  may  be  changed  at  the  execution  time  by: 


HNR52.  PUNCH=save  file. 

On  COC  computers  the  save  file  Is  TAPE?  by  default  and  may  be  changed 
at  execution  time  by: 


HHR52,  INflle,  OUTflle,  SAVEflle. 


On  MS-OOS  compatible  microcomputers  (PC's),  the  save  file  Is  UNIT8  by 
default  and  may  be  changed  at  execution  time  by: 


HMR52,  INflle,  OUTflle,  SAVEflle. 
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Hmi?  INPUl  OtSCSIPlION  . 

SA 

SA  CARD  -  STORH  AREA  AND  0R1EN1A110N 
FIEIJ  VARIABLE  VALUE  DESCRIPTION 

0  II)  SA  Card  Identification. 

1  PHSA  *■  Probable  maximum  storm  area  size  In  square  miles. 

If  PMSA  Is  zero,  the  program  will  compute  PMS  for 
several  area  sizes  and  select  the  area  size  which 
produces  maximum  precipitation  on  the  drainage 
basin  for  the  specified  number  of  6-hour  periods. 

2  ORNT  Probable  maximum  storm  orientation  In  degrees, 

clockwise  from  north.  (Range  Is  135-315  degrees.) 

If  ORNT  Is  less  than  or  equal  to  zero,  the  program 
will  compute  PMS  for  several  orientations  and 
select  the  orientation  which  produces  maximum 
precipitation  on  the  drainage  basin  for  the 

specified  number  of  6 -hour  periods. 

0  If  ORNT  Is  zero,  the  orientation  which  minimizes 
moment-of-lnertia  for  the  drainage  basin  about  the 
major  axis  will  be  used  to  estimate  PMSA  when  Field 
1  Is  zero. 

If  ORNT  Is  negative,  the  absolute  value  of  ORNT 
will  be  used  to  estimate  PMSA  when  Field  1  Is  zero. 

3  NINCS  f  Number  of  6-hour  periods  to  use  for  computing 

maximum  precipitation  on  the  drainage  basin 
(default  Is  3). 


II 


HHR5?  INPUT  DESCRIPTION 


SC 

gQ  SC  CARD  -  STORM  CENTER  CORRDINATES 

The  slorm  center  will  be  located  at  the  drainage  basin  centroid  If  an  SC 


card  Is 

not  used. 

FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

11) 

SC 

Card  Identification. 

1 

XCEN 

X-coordInate  of  storm  center. 

2 

YCEN 

f 

Y-coordInate  of  storm  center. 

SO  CARO  -  ARRANGEMENT  OF  6-HR  INCREMENTS  IN  PMS  TEMPORAL  DISTRIBUTION 

This  card  gives  the  arrangment  of  6-hr  Increments  of  precipitation  In  the 
PMS.  The  default  arrangement  (If  no  SO  card  Is  given)  Is:  12,  10,  8,  6,  4,  2, 
1,  3,  5,  7,  9,  11  where  the  numbers  Indicate  the  largest  (1)  to  smallest  (12) 
6-hr  Increment. 


If  SO  cards  are  Included  In  a  data  set,  all  twelve  Increments  must  be 
listed  on  two  cards . 


FIELD  VARIABLE 
0  II) 

1  RELMAG(l) 

2-10 


VALUE  DESCRIPTION 

SD  Card  Identification. 

f  Relative  magnitude  of  first  6-hr  period  of  PMS. 

*■  Relative  magnitude  of  remaining  6-hr  periods  of  PMS. 


Continue  with  eleventh  period  In  Field  1  of  second  SO  card. 
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HMR52  INPUT  DESCRIPTION 


SI  CARO  -  STORM  ISOHYET  AREAS 

This  card  is  optional.  If  SI  cards  are  used  the  program  will  not  calculate 
the  basin  area  within  each  isohyet,  but  it  will  use  the  areas  given  on  SI 
cards.  Storm  area  and  orientation  must  be  given  when  SI  cards  are  used. 

If  both  SI  and  BA  cards  are  used,  then  BX  and  BY  cards  may  be  omitted. 

If  SI  cards  are  used,  then  there  must  be  two  SI  cards. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

SI 

Card  identification. 

1 

AREAlB(l) 

Area  within  both  basin  boundary 
square  miles. 

and 

isohyet  A  in 

2-10 

AREAIB(2-10) 

f 

Area 

miles 

within  basin  boundary  and 
for  isohyets  B  through  J. 

isohyet  in  square 

SECOND 

SI  CARD 

riELD 

VARIABLE 

VALUE 

DESCRIPTION 

1-9 

AREAIB(11-19) 

¥ 

Area 

miles 

within  basin  boundary  and 
for  isohyets  K  through  S. 

isohyet  in  square 

10 

AREAIB(20) 

* 

Area 

in  square  miles  within  basin 

boundary  ^nd 

isohyet  coresponding  to  the  storm  area  size. 


SI 
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ST 


HNRS2  INPUl  DESCRIPTION 


ST  CARD  -  TENPORAL  DISTRIBUTION 


less 

This  card  gives 
than  6  hours. 

data 

for  the  temporal  distribution  of  the  PMS  for  Intervals 

FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

ST 

Card  Identification. 

1 

TIMINT 

f 

Time  Interval  In  minutes  to  be  used  for  temporal 
distribution  of  PMS.  (Range  of  5  min  to  360  min). 

2 

R16A20 

1- 

Ratio  of  ihr  to  6  hr  precipitation  for  Isohyet  A  of 
20,000  sq  ml  storm  from  Figure  39  of  HMR  No.  52. 
(Range  0.2?  to  0.35). 

3 

POSMAX 

Position  of  maximum  6-hr  Increment  In  PMS  temporal 

distribution.  Remaining  6-hr  Increments  In 
descending  order  will  be  placed  alternately,  before 
and  after  maximum  Increment.  This  arrangement  will 
replace  any  distributions  from  previous  SD  or  ST 
cards.  (Range  5  to  12,  default-7). 

0  Use  previously  established  arrangement  of  6-hr 
Increments  of  PHS. 

4  RATIO  ♦  PMS  precipitation  will  be  multiplied  by  this 

ratio.  (Default  =  1.0.) 
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HMR52  INPUl  DESCRIPIION 


zw 


ZW  CARD  -  WRITE  PMS  TO  OSS 

The  ZW  card  Is  IncTuded  with  data  for  each  basin  for  which  the  hyetograph 
Is  to  be  saved  on  OSS. 

Each  data  record  (hyetograph)  Is  Identified  by  a  6-part  pathname,  each  part 
being  designated  by  Tetters  A  through  F. 

Part  A  Is  project  Identification  or  description,  T6  characters  maximum. 
Default  Is  blank. 

Past  B  Is  hyetograph  location,  8  characters  maximum.  Default  1s  basin  name 
from  BN  card.  Field  T . 

Part  C  Is  parameter  name,  12  characters  maximum.  Default  Is  PRECIP-INC. 

Part  D  Is  start  date  of  hyetograph.  Default  Is  01JAN1999.  Hyetographs  are 
stored  In  blocks  of  1  month  or  1  day  depending  on  time  Interval.  The  hyetograph 
will  start  at  0000  hours  on  the  date  given  on  the  ZW  card,  but  the  date  In  the 
pathname  will  be  the  beginning  date  of  the  block. 

Part  E  Is  time  Interval.  This  will  be  derived  by  the  program. 

Part  F  Is  alternative  or  description  of  hyetograph,  24  characters  maximum. 
Default  Is  blank. 

The  format  of  the  ZW  card  Is 

ZW  A=Part  A,  B=Part  B,  C=Part  C,  0=Part  D,  F=Part  F 

One  or  more  of  the  pathname  parts  may  be  omitted  from  the  ZW  card.  Once 
parts  A,  C,  0,  or  F  have  been  set  they  will  retain  their  values  until  reset  on  a 
ZW  card. 

For  example: 

ZW  A  =  F=PNS-1 

Will  set  part  A  to  blank,  and  part  F  to  PMS-1. 
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HHR5?  INPUT  DESCRIPTION 

a  CARD  -  END  OF  DATA 

FIELD  VARIABLE  VALUE  DESCRIPTION 

0  II)  //  Card  Identification. 


HMR52  INPUT  DESCRIPTION 


EXAMPLE 

II)  PROBABLE  MAXIMUM  STORM  CALCULATION  FOR  EXAMPLE  2  IN  HMR  NO.  52 

ID  OUACHITA  RIVER  BASIN 

*fR£E 

BN  OUACHITA 

**  CALCULATE  STORM  OVER  ENTIRE  BASIN 
BS  .79365 

BX  108  108  107  97  98  97  95  93  88  83  77  71  67  66  64  62 
BX  59  54  48  40  38  39  37  26  18  12  8  6  4  6  8  9  15  25  27 

BX  30  33  35  38  41  47  53  56  67  70  76  80  82  85  90  100  105 

BY  17  21  24  32  34  36  37  44  48  48  46  46  50  54  55  54  50  47 

BY  47  43  41  38  37  42  39  40  40  42  41  37  34  28  23  19  20 

BY  20  19  23  21  21  24  22  23  16  15  10  10  14  11  12  10  12 

HP  10  30.0  35.9  40.6  44.6  47.1 

HP  200  22.2  27.0  31.2  34.7  37.7 

HP  1000  16.3  21.0  25.3  29.0  31.2 

HP  5000  9.5  13.5  17.7  21.6  24.2 

HP  10000  7.3  10.7  14.0  18.0  20.0 

HP  20000  5.3  8.5  11.6  14.9  17.2 

HO  235 
** 

**  FIND  STORM  AREA  AND  ORIENTATION  FOR  PMS  CENTERED  AT  BASIN  CENTROID 
SA  0  0 
GO 

**  CALCULATE  STORM  CENTERED  BEIWEEN  RENNEL  DAM  AND  BLAKELY  MT.  DAM 
**  USING  GIVEN  STORM  AREA  AND  ORIENTATION 
SC  87.9  20.1 
SA  2150  280 
GO 

**  CALCULATE  S10RM  FOR  EACH  SUBBASIN  USING  LAST  SEl  OF  PARAMETERS 
BN  PINERIDG 

BX  26  18  12  8  6  4  6  8  9  15  25  27  30  33  35  33  35  33  35 
BY  42  39  40  40  42  41  37  34  28  23  19  20  20  19  23  27  31  35  38 
BN  NASHUA 

BX  35  38  40  40  44  47  57  63  58  55  54  51  53  48  40  38  39  37  35  33 
BX  35  33 

BY  23  21  21  23  25  30  31  35  42  42  44  44  47  47  43  41  38  37  38  35 
BY  31  27 
BN  BLAKELY 

BX  53  51  54  55  58  63  57  47  44  40  40  41  47  53  56  59  64  72  79  84  85 

BX  96  95  93  88  83  77  71  67  66  64  62  59  54 

BY  47  44  44  42  42  35  31  30  25  23  21  21  24  22  23  21  23  22  28  28  30 

BY  36  37  44  48  48  46  46  50  54  55  54  50  47 

BN  RENNEL 

BX  108  108  107  97  98  97  96  85  84  79  72  64  59  67  70  76  80  82  85  90  100 

BY  17  21  24  32  34  36  36  30  28  28  22  23  21  16  15  10  10  14  11  12  10  12 

ZZ 
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APPENDIX  B 


Error 

Number 

1 


2 


3 


ERROR  MESSAGES 


Error  Message 
and  Description 

STORM  AREA  AND  ORIENTATION  MUST  BE  GIVEN  WHEN 
ISHOYET  AREAS  ARE  GIVEN  ON  SI  CARDS 

Storm  area  Is  required  to  Identify  the  division  between 
wlthln-storm  and  residual  rainfall.  Storm  orientation  Is 
required  to  calculate  orientation  adjustment  of  PMP. 


BOTH  BA  AND  SI  CAROS  MUST  BE  USED  WHEN  BOUNDARY 
COORDINATES  ARE  NOT  GIVEN. 

The  program  cannot  calculate  basin  area  or  areas  enclosed 
by  Isohyets  and  basin  boundary  without  boundary  coordinates. 
These  values  must  be  given  on  BA  and  SI  cards  If  boundary 
coordinates  are  not  given. 


AREA  WITHIN  BASIN  DOES  NOT  INCREASE  WIIH 
ISOHYET  SIZE.  AT  =  xxxx.xx.  A?  =  xxxx.xx 

This  error  occurs  when  the  program  Is  computing 
basin-average  precipitation  and  the  area  within 
the  basin  and  an  Isohyet  Is  smaller  than  the  basin 
area  Inside  the  next  smaller  Isohyet. 

This  error  may  be  caused  by  an  Incorrect  entry  for 
the  basin  area  Inside  the  storm-area  Isohyet  on 
an  SI  card. 

This  error  may  also  be  caused  by  Incorrect  boundary 
coordinates.  If  the  basin  boundary  crosses  Itself 
the  program  may  compute  a  wrong  basin  area  Inside 
an  isohyet. 


BX  CARO  NO.  nnn 

8X  CAROS  HAVE  ALREADY  BEEN  READ  FOR  THIS  BASIN. 

A  BN  CARO  MAY  BE  HISSING,  OR  A  CARO  MAY  BE  OUT  Of  ORDER. 

BX  cards  for  a  single  basin  may  not  be  Intermingled 
with  other  card  types.  A  BN  card  Is  used  to  Identify 
the  beginning  of  data  for  a  basin. 


8t. 


BY  CARD  NO.  nnn 

BY  CARDS  HAVE  ALREADY  BEEN  READ  FOR  THIS  BASIN. 

A  BN  CARD  MAY  BE  MISSING.  OR  A  CARD  MAY  BE  OUT  OF  ORDER. 


BY  cards  for  a  single  basin  tnay  not  be  Intermingled 
with  other  card  types.  A  BN  card  Is  used  to  Identify 
the  beginning  of  data  for  a  basin. 


BA  CARO  NO.  nnn 

BASIN  AREA  ON  BA  CARD  IS  NOl  GREATER  THAN  ZERO. 

AREA  =  xxxx.xxx 

Basin  area  must  be  greater  than  zero.  If  program  Is  to 
calculate  area  from  boundary  coordinates  BA  card  should  not 
be  used. 


HP  CARD  NO.  nnn 

DEPTH  DOES  NOT  INCREASE  WITH  DURATION  FOR  AREA  =  xxxxx. 

PMP  depths  on  the  HP  card  for  area  xxxxx 
should  Increase  from  field  ?  to  6. 


INVALID  AREA  ON  HP  CARD  NO.  nnn,  PMP  AREA  =  xxxxx. 

PMP  area,  field  1  of  HP  card,  must  be  10,  200,  1000, 
5000,  10000,  or  20000  square  miles. 


SO  CARD  NO.  nnn 

MAGNIIUDE  (m)  IS  NOl  IN  RANGE  1  TO  12. 

Values  on  the  SO  card  are  Integers  from  1  to  12 
which  Indicate  the  relative  magnitude  of  6-hour 
precipitation  Increments.  "1"  Indicates  the  largest 
Increment  and  "12"  Indicates  the  smallest  Increment. 


SD  CARO  NO.  nnn 

A  RELATIVE  MAGNIIUDE  WAS  REPEATED 
Ihe  values  on  the  SO  card  are  unique  Integers 
from  1  to  12.  No  value  may  be  repeated  on  a 
set  of  SD  cards . 


SD  CARD  NO.  nnn 

TWO  SD  CARDS  ARE  REQUIRED. 

Twelve  values  are  read  from  SO  cards.  This 
requires  two  SD  cards  In  sequence. 


SI  CARD  NO.  nnn 

AREA  FOR  ISOHYEl  11  (nl)  IS  LESS  THAN  AREA 
FOR  ISOHYEl  11 

Areas  enclosed  by  Isohyets  and  the  basin  boundary 
cannot  decrease  as  Isohyet  area  Increases. 


ST  CARD  NO.  nnn 

TIME  INTERVAL  CANN07  BE  LESS  THAN  5  MINUTES 

The  program  will  not  calculate  a  temporal  distribution  for 
a  time  Interval  smaller  than  5  minutes. 


SI  CARD  NO.  nnn 

TIME  INTERVAL  MUST  DIVIDE  60  MINUTES  INTO  EQUAL 
PARTS  WITH  NO  REMAINDER 

The  time  Interval  for  temporal  distribution  must 
be  5,  6,  10,  12,  15,  20,  or  30  minutes. 


ST  CARO  NO.  nnn 

TIME  INTERVAL  HAY  NOT  BE  GREATER  THAN  6  HOURS  (360  MINUTES) 

The  program  will  not  calculate  a  temporal  distribution  for 
a  time  Interval  greater  than  6  hours. 


ST  CARO  NO.  nnn 

TIME  INTERVAL  MUST  DIVIDE  6  HOURS  INIO  EQUAL 
PARTS  WITH  NO  REMAINDER 

The  time  Interval  for  temporal  distribution  roust 
be  1 ,  2,  3,  or  6  hours . 


ST  CARD  NO.  nnn 

1-HR  TO  6-HR  RATIO  FOR  ISOHYEl  A  OF  20000  SQ.  MI. 
STORM  MUSI  BE  IN  THE  RANGE  0.27  10  0.35 

The  range  of  the  1-hour  to  6-hour  ratio  from 
figure  39  of  Hydrometeorological  Report  No.  52 
Is  0.270  to  0.300.  A  value  outside  this  range 
Is  considered  to  be  a  mistake. 


SI  CARO  NO.  nnn 

POSH  ION  Of  MAXIMUM  6-HR  INCREMENT  MUSI 
BE  IN  RANGE  5  TO  12 

The  maximum  6-hour  precipitation  Increment 
may  not  occur  In  the  first  24  hours  (positions 
1  to  4)  of  a  storm. 


19  ARRAY  ON  CARD  NO.  nnn 

(id)  EXCEEDS  DIMENSION  OF  kkk 

The  program  attempted  to  read  more  data  than 
allowed  from  an  xx  card.  The  maximum  number  of 
values  that  can  be  read  from  xx  cards  Is  kkk. 


n  ERRORS  DETECTED  IN  INPUT  DATA  FOR  BASIN  aaaaaaaa. 
REMAINING  DATA  WILL  BE  CHECKED  FOR  ERRORS. 

n  ERRORS  DETECTED  IN  CALCULATIONS  FOR  BASIN  aaaaaaaa. 
REMAINING  DATA  WILL  BE  CHECKED  FOR  ERRORS. 

When  an  error  Is  detected  by  HMR52,  the  program 
will  read  through  the  remaining  data  and  check 
for  Input  errors.  No  calculations  will  be  made 
using  the  remaining  data. 
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Computer  requirements: 
lime  and  memory,  6? 

File  characteristics,  64 
Error  messages,  85 
Flood  (see  Probable  Maximum  Flood) 

Hydrometeorological  Report  No.  52  (HMR  No.  5?) 

Nomograms  for  PMP  Increments,  25 
Input  data: 

Description,  67 
Overview,  40 

Output  from  HMR52  program: 

General,  42 

Precipitation  file,  49 

Precipitation  (see  Probable  Maximum  Precipitation  or  Storm) 
Probable  Maximum  Flood  (PMF): 

Determination,  29 
Example  calculation,  54 
Parameter  sensitivity,  55 
Probable  Maximum  Precipitation  (PMP): 

Definition,  3 
Determination,  14 
Example  calculation,  15 
Regional  guidance,  1 
Probable  Maximum  Storm  (PMS): 

Basin-average  precipitation,  35 

Definition,  3 

Examples: 

Hand  calculation,  18 
Computer  calculation,  50 
Optimize  storm  area  and  orientation,  39 
Orientation,  4,  33,  34 
Shape,  3 
Size,  4,  13,  31 
Spatial  variability,  6,  9 
lemporal  distribution,  7,  36 
Watershed  and  storm  geometry,  30 
References,  65 

Runoff  (see  Probable  Maximum  Flood) 

Watershed  geometry,  30,  34 
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